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ABSTRACT 

Vertical  stern  accelerations  for  five  geometri- 
cally similar  destroyer  hull  forms  of  different  lengths 
in  the  250-550  foot  range  are  predicted  on  the  basis  of 
model  tests  for  a  ship  speed  of  17  knots  at  various 
headings  in  representative  short-crested  irregular  seas 
of  varying  severity. 

Using  trie    concept    of   statistical   inference,    the 
model  test   predictions   are    treated   as   a   statistical 
sample    and   used   to   predict   what  would   occur   in   the    total 
population  of   possible    stern  accelerations. 

Then,    with  strong  emphasis    on   the    systems 
approach,    it    is    shown  how   these    predictions   can   be    used 
as    input   to   a   probability   analysis    that    leads    to   a   pre- 
diction  of   long  range    system  behavior/effectiveness, 
the    system  under   consideration   being  the   Destroyer/ 
Variable   Depth  Sonar  System. 


DD  -  Destroyer  type  ship. 

VDS         -  Variable  Depth  Sonar. 

RAO         -  Resoonse  Amplitude  Operator. 

y)  -  Actual  wave  length  of  regular  waves  in 

feet  measured  perpendicular  to  lines  of 
crests . 

i 

Effective^     -  Wave    length  apparent    to   a   ship   proceeding 
obliquely   through  regular  waves. 

L  -   Ship   length   in   feet    measured   along 

designer's   waterline 

Amplitude        -  For  a  regular   or    irregular  function, 
peak  to    peak  value    civided  by   two. 
Sometimes   referred   to  as    single 
amplitude • 

Height  -  For   a  regular    or   irregular  function, 

peak  to   peak  value.      Some  tires   referred 
to   as   double    amplitude. 

u>  -  Amplitude    of   ship   pitch   angle    in  radians; 

'*  vertex   of   angle    is   defined   as   being   at 

ship's   center   of   gravity. 

^  -   Instantaneous   value    of    pitch. 


<c 


Z0  -  Amplitude  of  heave  (vertical  transla- 

tion) of  shin  in  feet  at  center  of 
gravity. 

Z0  -  Instantaneous  value  of  heave. 

Si  -   Amidships 

2X  -  Vertical  translation  at  a  point  x  feet 

from  Q. 

h  -  Wave  amplitude  in  feet. 


ii 


h  -  Instantaneous  value  of  distance  In  feet 

between  wave  surface  and  calm  water  zero 
datum* 

\}max       -  Maximum  wave  slope.   See  Page  9* 

^_  -  Phase  angle  between  pitch  and  wave. 

~V  See  Pages  8  and  9. 

-  Phase  angle  between  heave  and  wave* 
See  Pages  8  and  9* 

£         •  Phase  angle  between  heave  and  pitch* 


it- 


See    Page   9* 


ycy    „  /i  -  See  Figure  3  on  page    11  for  angular 

'   .S*,/^*»>  nomenclature. 

v8  -  Ship  speed   in  feet   per  second. 

l>j  -  Circular  frequency  of  wave    in  radians   per 

second* 

u^  -  Encounter  frequency  in  radians   per  second* 

t  -  Time    in  seconds. 

g  -  Acceleration  due   to  gravity;  32.2  ft.   per 

(second )« 

x  -  Distance  along  ship's    length  from  8  to 

some   point,    (+)  forward,    (-)  aft.      In 
Appendix  A,   ordinate    of   irregular  stern 
acceleration  record. 

Ax  -  Amplitude   of  vertical  acceleration  at  a 

point  x  along  ship's    length.      In   this 
thesis  x  =  -181  feet  for  the  370  foot 
ship  length,      (x  changes  with  geometric 
similarity  for  other  ship  lengths) 

In  -  Logarithm  to  the  natural  base  e. 

(T  -  Increment   of. 

corr.ft*  -   Pierson-Moskowitz   sea   spectrum  ordinate. 

"Corr"    implies  corrected.     See   Page    16 

H  -  Lag  number;   sea  spectrum  abscissa  fre- 

quency coordinate.      See    Pages   16  and   17. 


ill 


Hw3  -  Significant  wave   height    in  feet;  defined 

as   the   average    of   the    1/3  highest  wave 
heights  extant   in  any  given  seaway. 

RSC  -  Response   Spectrum  Component.     See  Figure   9. 

IRS  -  Integrated  Response  Spectrum.     See  Figure   9. 

l/E  -  Ship  response    parameter.      See  Appendix  A. 

Sample  -  A   relatively  small  portion  of  the    set   of 

values   in  the   universe   of  a  random  variable. 

x  -  Mean  of   sample   probability  distribution  of 

stern  vertical  ^JE  acceleration. 

Q 

8  -  Variance   of   sample   probability  distribution 

of   stern  vertical  tE  acceleration. 

s  -  Deviation  of   sample   probability  distribu- 

tion of  stern  vertical  "tfE  acceleration. 

N(x,s    )  •  A  normal   (Gaussian)  probability  distribu- 

tion function  whose  mean   is  x  and  whose 
variance   is  S2. 

Population     -  The   set  of  all  possible  values   in  the 
universe   of  a  random  variable. 

f*  -  Mean  of  population  probability  distribution 

of  stern  vertical  "\JE  acceleration. 

<T   2  -  Variance   of  population  probability  distri- 

bution   of  stern  vertical  liE  acceleration; 
when  used   in  a  general  sense,    the  variance 
of  any  normal  probability  distribution 
function. 

f  "  Deviation  of  population  probability  distri- 

bution of  stern  vertical  "tfE  acceleration* 

^o^o  -  True   population  parameters. 

Pr  [    1  -  Probability  of  [  ] . 

Au,  (f\x  -  Estimate   of  ^b,  ^o.     See  Figure    17  on   Page  56. 

P  -  Location  of  x  and  8   in  <J^v  parameter  space. 

n  -  Size   of  a  statistical  sample. 


iv 


N  -  Number  of  cycles  of  a  certain  ship  motion. 

R  -  Parameter  of  the  Rayleigh  probability  dis- 

tribution function. 

RMS         -  Root  mean  square  of  some  variable 0   See 
Appendix  A. 

1/fT        -  See  Appendix  A. 

Superscripts 

-  Denotes  amplitude  of  a  sinusoidal 
function  as  distinguished  from  instan- 
taneous value.   In  the  case  of  an  irregular 
function,  denotes  peak  to  peak  divided  by 
two.   This  superscript  does  not  imply  the 
above  when  used  in  x;  in  that  case  it 
indicates  an  arithmetic  mean. 

-  Average    of. 

Subscripts 

9  -  Measured   at   ship's   center   of   gravity. 

max  -  Maximum  value. 

2*  -  heave   relative    to  v/ave  . 

y  -   pitch  relative    to  wave. 

ZV  -   pitch  relative    to  heave. 

x  -  at    location  x  feet   along  ships    length 

from  B. 

e  -  Encounter. 

^  -  the  i  th  value  of  a  variable. 

u  -  upper. 

L  -  lower. 

PTP  -   peak-to-peak  value    (double    amplitude    or 

height ) 


l/5       -  the  highest  l/3  of  a  group  of  values 
of  a  variable . 

1  -  original  value,  before  being  changed, 

2  -  value  after  having  been  changed  from 

original  value. 
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INTRODUCTION 

The   basic  method   used   in  this   thesis   to  predict 
the  various  modes   of  response  of  a  particular  ship  to  a 
given  seaway  was   introduced  over  ten  years  ago    (1)#  and 
has  been  extensively  employed   in  the    intervening  time* 
The   challenge   of  the   present  day  is  how  to  most  effec- 
tively utilize   such  response  data   in  order  that   the 
systems  de s igne r  might  be   provided  with  reliable    long- 
range  estimates   of  system  behavior  and  as  well,   that   the 
system  user  might  have   reliable  estimates   of  system 
effectiveness* 

The   approach  developed   in  this   thesis   is  a 
broad  and  general  one   that    is   applied  to  one    particular 
aspect   of   behavior  of   one   system*      It   is   suggested  that 
the   reader  keep  in  mind  the   generality  of  the  method, 
while   considering  the   specific   case  which  is  developed 
herein  for  illustration  as  well  as  for  its   current 
naval  interest* 

The   system  considered   in  this  thesis   is   the 
DD/VDS  system.     The   particular  aspect  of  system  behavior 
considered   is  the  vertical  acceleration  of  the  YDS  tow- 
point  at  the  destroyer^   stern* 


^Numbers    in  parentheses  refer  to  references   listed   on 
Page  91. 


The   development    undertaken  herein  can   be   divided 
into  three   principal  areas,  namely! 

(1)  prediction  of  vertical  stern  acceleration 
(in   terms   of   the  'vS  parameter  which  is 
discussed   in  detail  in  Appendix  A),  for 
various   ship  lengths,  headings  and   sea 
states; 

(2)  utilization  of  the   concept   of  statistical 
inference   to  predict   the   population  dis- 
tribution of  vertical  stern  acceleration* 
(in  terms   of  the  ^  E  parameter),   treating 
the   predicted  values   in    (1)   above  as  a 
statistical  sample; 

(3)  a  demonstration  of  how  these   predicted  stern 
accelerations   can  be    used   in  a  probability- 
based   system  analysis    leading  to   long- 
range   prediction  of  CD/VDS  system  behavior 
over  the    lifetime    of  such  a  system* 


*The   term  "population"   is   used  here   to  signify  all  the 
possible  values   of  vertical  stern  acceleration    (in  terms 
of  the  V!  parameter)   that  a  destroyer  of  the    particular 
type  considered  might  experience   in  its   lifetime.      In 
theoretical  probability  and   statistics,  "population"   and 
"universe"   are   used   interchangeably  to  describe   the   set 
of  all  the   possible   values   of  a  given  random  variable. 


PREDICTION   OF   STFRN  ACCFLF1- ATTON  RES PONSE 

General ^Procedure    for    Prediction    of    Ship  Response 

In   1953,    St.   Denis    and    Pierson    (1)    introduced    a 
method   utilizing  the    principle    of   superposition  whereby 
a   ship's   response    in   irregular  waves   may   be    preoicted 
provided    the    response    in  regular  waves    is    known,    as    it 
is,    for  example,    from  model   tests    in   regular  waves. 

The   method   basically  consists    of    (a)   representing 
the    sea   as    an  energy   spectrum,    (b)    transforming  model 
response   data    taken    in  regular  waves    into  a  response 
amplitude    operator    (RAO)    (referred    to   by   some    as   a 
transfer  function)   and    (c  )   multiplying  the   RAO  by   the 
sea   spectrum  to  arrive    at   a   response    spectrum  having 
certain   probability   properties   regarding  the   response 
being  considered. 

One    section  each  will  be    devoted    to   the   descrip- 
tion  of  developing  RAO's   from  model  data,    to   the    selection 
of   sea   spectra,    and    to   the    computational  procedure    used 
to  combine    sea   spectrum  and   RAO   to   arrive    at    a  response 
spectrum. 

Development   of  RAO's   from  Model  Data 

The   most   comprehensive    source    of  destroyer  model 
motions   data  available   at    the    present   time    is   the  Dutch 


data   that   can   be    found    in    (2).      It    is    particularly- 
unique    in   that   model  response    at   various  headings    to 
the   waves   was   measured   and    this    is,    of   course,    a 
necessary   prerequisite    if   one   desires    to  make    predic- 
tions   in   short-crested    seas   and   to   investigate    the 
general  effect    of  heading   on   ship  response.      This  Dutch 
data   is    also   unique    in    that    phase    angles   were    measured 
and    presented.      The    primary   limitation   of    this   data   is 
that    the    above    variations   were   done    at    only   one    speed, 
namely   17   knots    full   scale. 

The    vertical   acceleration   at   the    stern   of   a   ship 
is   caused   primarily  by   pitching  and    somewhat    less   by 
heaving.      Roll   is    the    only   other    possible    contributor 
to  vertical  acceleration,    but    its   effect   was    ignored 
because    the    stern   VTS   towpoint    is    on   the    shipfs   center- 
line.      The    vertical   acceleration   also    is   a   function   of 
the    phase    angle    between  heaving   and    pitching.      Conse- 
quently,   the    data  from    (2)    used    in   this    investigation 
consisted    of  heave,    pitch  and    the   heave-to-wave    and 
pitch-to-wave    phase   angles.      Samples    of   the   model  data 
are    shown   in   Figures    1  and   2. 

The  model  used  was  17002.  It  is  a  model  of 
Her  Netherlands  Majesty's  Ship  Priesland.  Its  hull 
parameters   are    as    shown   in  Table    I. 
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TABLE   I 


L  370  ft.  Waterplane   Coefficient    0.806 

Beam  33.6  ft.  Prismatic   Coefficient      0.688 

Draft  13.2  ft.  Displacement  3,010  tons 

Block  Coefficient  0.566  Number  of  Rudders  2* 


It   can   be    seen  from  Figures   1  and  2   that  a  certain 
amount  of  extrapolation  of   the   model  data  was  necessary. 
On  the   high  end   of  the  x/L  scale,   2Q->h  and  Y0-*2Tl  hA  s 
Uhnax,    but   the   question   immediately  arose   as   to    just 
where   on  the    Vl  scale   it  would   be   reasonable   to  assume 
that    this   occurs.     Extensive  examination  of   other  destroyer 
model  data    (3),    (4)   showed   that  a  */L  of  3»75   to  4-00  is 
reasonable  for  heave   and   pitch  having  attained  their 
limiting  values.      Since   the   angle    between   ships   course 
and  wave  direction  had   to  be   taken   into  account,  T0  was 
faired   to  effective  wave   slope    (based   on  effective  A, 
where  effective  A  is  the  A apparent  to  a   ship  proceeding 
obliquely  through  the  waves)  and  2Q  was  faired   to  h  at 
whatever  actual  *  resulted    (about  3.00  in  the  cases   shown 
in  Figures  2  and  3)   when  effective  */L-3.75  to  4.0O. 

As  for   phase  angle  extrapolation,    the   statements 
in    (2)  are  generally  supported  by  vibration  theory  for 
a  single   degree   of   freedom  system.     The    limiting  values 
to  which  £r  and   jL4were  extrapolated    (using   the    same 
effective  7L limits   as  described  above)  were   taken  from 


the  following: 

£         ranges  from   (-180°  to  0°) in  head, 
fc       bow  and   beam  seas  as  effective   A/L 
increases   from  zero  to  about   4.00, 

£        ranges  from   (0°  to  0°)   in  quartering 
h       seas   over  the   effective   A/L  range   of 
zero  to  4.00. 

£         was   zero  throughout   the   T/L  range  for 
*       following  seas   because  of  very  small 
encounter  frequency. 

£         ranges  from  (+90°  to  -  90°)   in  head  and 
bow  seas  as  effective  7\/h  increases  from 
zero  to  about   4,00. 

£       was  90°  throughout   the    tyL  range    in 
beam  and  following  seas.      In   the   beam 
seas,  A  »      ship  beam  and  thus   there 
is   effectively  zero  pitch  phase    lag 
and  in  the  following  seas   case  ue    is 
very  small* 

£        ranges  from  (90°  to  90°)   in  quartering 
seas. 

Before   proceeding  further,   it   is   in  order  to  pre- 
sent  the   basic   conventions  and   relationships   used   so  as 
to  prevent   confusion.     There   are  many  different    systems 
in  use   so  the  best   policy  obviously  was  to  set  up  a 
system  and  adhere   to  it.      I  have  followed  that   policy 
based  on  these   convent icns  and  relationships:* 

(a)     If  the   wave   motion  amidships   is 
described   by 

h  ■  h  cos    (>*Je  t ) 

the  resulting  pitch  motion  is  described  by 

Yo  -  ^o  cos(^e  t  +£Y) 

♦In  the  notation  used  here  a  bar  (   )  above  a  symbol  im- 
plies an  amplitude  as  distinguished  from  an  instantaneous 
value. 
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and   the  heave  motion  by 

ZQ  a  Z0  oos(tJe   t   ♦  £*) 

(b)     heave   is   taken  as   positive   upwards 
and 
positive   pitch  angle   is   taken  as  bow 
downwards* 

(o)     it  can  be   seen  from  the  equations   in 
(a)   that  negative   phase  angles   imply 
motion  lagging  wave  and   positive 
phase  angles   imply  motion  leading  wave 

(d)  The   phase  angle   between  heave  and  pitch 
is  defined  by  the  equation 

Negative  values   of  £ir  imply  heave 
lagging  pitch  and  positive  values 
of  £      imply  heave    leading  pitch. 

(e )  Prom  basic  wave   theory  we   know  that 
maximum  wave   slope 

lLax  =££k 

(f )  It  can  be   shown  that  encounter 
frequenoy 

i  q 

g 


where     ^a  f2"TTg/^.  according  to  basio 
wave  theory. 

(g)     The  Butch  models  were   tested   in  waves 
of  proportions : 

Equipped  with  (a)  through  (g)  above,  it  is 
possible  to  determine  the  amplitude  of  vertical  dis- 
placement at  any 


*See  Figure  3  for  angular  nomenclature. 
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point  x  along   the    ship's    length,    Zx,    from 

Zx2    =  Z02    +  x2y02    .  2I0x  rocos   £*r  (5) 

where   x  =  distance    along   snip's    length  from  amidships    to  the 
point   under   consideration,    (+)   forward    (-)   aft.      In   this 
thesis   x  =  -181  feet   for   the   370  foot    ship  length  and  x 
changes    in  a   geometrically   similar  manner  for   the    other 
ship   lengths    investigated.      The   x  distance   was    chosen 
to  be    a  reasonable    location   of   the   VDS   towpoint. 

The    amplitude    of  vertical  acceleration  at   any   point 
x  along  the    ship's    length,   Ax,    can    be    computed  from 

Ax  =  WeS   Zx 
The    ordinate    of   the   RAO   is 

2  it  HA  J 
and    tne    abscissa   is    In  w. 

RAO's   for^    -   180°,    135°,    90°,    45°and    0°  were 
computed  directly  from   the   model  data   curves,   extrapolated 
as   necessary   to  cover   a   In  W  range    of   +0.2   to   -1.6.      The 
reason  for   the    limits    of   this   range   will  become   evident 
later. 

Then,    in   order   to  have   RAO's   at  22.5°   increments, 
the    original  RAO's   were   first   plotted   on  an  OJfc   base    to 
align  the    peaks   and    then  crossplotted   at   constant   values 
of  6)e.      The   resulting  RAO's   f  or  f4   -   157   l/2°,    112    l/2°, 
67   1/2      and   22    l/2      were    then  converted    to  a   In  dJ  base 
and    plotted   with   the    original  ones,    resulting   in  a 
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complete    selection   of  RAO's   at   22    l/2°   Increments   as 
shown   in  Figure   4. 

The   RAO's    thus   derived  are   for  the    full   scale 
length   of  the    destroyer  hull  considered,    namely  370  feet 
(this   was    the    L  used    in   the  3./L  scale    of   the    original 
model  data ) . 

To   arrive    at  RAO's   for  geometrically   similar 
hulls    of   other    lengths    is    a   simple   matter  when   they  are 
plotted    on  a   In  W  base    (5).        As    stated   previously, 


W   ="|/2?g/A 
At  7)/L  a  a   constant    value, 


U^  =  k-jjl/Li  Subscript    1   indicates 

original  value 


Wg   =  k  IJI/L2  Subscript   2   indicates 

changed   value 


and 


ui  = 


4  V  wi-i 

Taking  the   natural  base    logarithm  of  both   sides 
of   this   equation,    In  U)^   -  In  Ug   =  l/2(ln  Ig    -   In  Li )    and 
therefore,    In  W2   -   In  Wi  =    6  In  OJ  =  -1/2 (In   L2   -   In   Li). 
It  follows  that:      c^ln  fc  =  -l/2   lnj1^} 

This    last  equation   shows    that  RAO's   for   other 
ship   lengths   can  be    obtained   by   translating  the    original 
RAO's   along   the   abscissa. 
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For  example,    let    cfln  u  ■  0.1 


0.1   =   -1/2    In   Lg 

*4      =    .819 
Ll 

and    Lg    =    ,819    (370)    -  502' 

Similarly,    if  (fin  10  =  -0.1,      Lo 

r=-  =  1.222   and   Lo   =  452 ■ 

Thus,    if   the   RAO ' s   are    shifted   to   the   right,    keeping  the 
In  &  scale    fixed    in    space,    ship   length   increases    and 
conversely  a   leftward    shift   causes   a  decrease    in  ship 
length.      Shifting   increments    of   0.2,    0.1,    -0.1,    and 
-0.2   were    chosen,    resulting   in  RAO's   for   snip   lengths 
of   248  ft.,    5  02    ft.,    452    ft.    and    552   ft. 

Figure    5   shows    the    157   l/2°  RAO   shifted   as 
described   by  the   above    procedure.      All   the    others   were 
shifted    in  exactly    the    same   manner. 

Selection   of   Sea   Spectra 

Through  the    principle    of   superposition,    the 
large   number   of  component  waves   comprising  an   irregular 
long-crested    seaway  may   be   represented   in  a   spectral 
form  having   an   ordinate   representing  energy  and    an 
abscissa  representing  wavelength   or  frequency.      This    is 
a  well-established   concept   and    is   described  rather  well 
in    (6).    (The   term  long-crested   is   used   by  oceanographers 
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to  describe    the    situation  where   all  the  waves   in  a  system 
are   coining  toward   the    observer  from  the    same   direction 
such  that  he    sees   a  well-defined   pattern  of   long,    parallel 
crests  )• 

The    log-slope   form  of   sea   spectrum  representa- 
tion developed  by  Lewis    (5)  was   used   in   this    present 
work.      For  my  purpose    it  has    one    large    advantage    over 
other  representations   which  have,    for  example,    (wave 
amplitude )2      as    ordinate   and  encounter  frequency  or 
lag  number  as   abscissa,   namely: 

RAO's   for  geometrically  similar  hulls 
of  different    lengths   can  be   easily 
developed   by  the    simple    expedient   of 
translating   them  along  the    In w>  axis* 


rUnuJ 


The   ordinate    of  the    log-slope   sea   spectrum  is 
2 


and    the    abscissa   is    lnw, 


Since   a   statistical  sampling  approach  was   de- 
sired  for   this   analysis,    the    Pierson-Moskowitz   collec- 
tion of   sea   spectra    (7),    (8)   was   chosen  because    it    is 
the    largest,   most   recent   collection   of  actually  measured 
sea  spectra.      Pierson  and  Moskowitz  use   "corr.ft2"*  as 
their  spectral   ordinate    and   lag  number,   H,   for   their 
abscissa.      See   Figure   5A  for  an   illustration   of  a   typical 
Pierson-Moskowitz   sea  spectrum.    The    pertinent   conversion 


*Corr,ft2  means  "corrected"   ft2-Pierson  and   Moskowitz  cor- 
rected  some    spectral  ordinates   to   smooth  out    occasional 
unreasonable   values   due   mainly   to  vagaries   in  the   measured 
waveheight  data. 
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relationships   are    as   follows: 

H  =   180  x   frequency   =180  ^ 

2  n 


r(ln  uj) 
*/2  TT~ 


I2   c   ^ 


r(  ui) 


=  i^£     x  180  x    (corr.ft2) 


A   short   computer   program,    written   in  Fortran    II 
source    language,    and    shown   in  Figure    6,    was    used    to   con- 
vert   the . Pierson-Moskowitz   Spectra   to   the    log-slope    form 
used    in    this    bhesis.      A   sample    of    the    print -out    is 
shown    in  Figure    7. 

The    next   matter   of    interest    is    that    of   short- 
crestedness.      This    is   not    an  easy  matter    to  dispose    of, 
simply  because    a   short-crested    sea  may  arise   from 
different   causes.      As   its  name    implies,    a   short-crested 
sea   is   one   where    the   waves  are   coming  toward   the    obser- 
ver from  various   directions,    thus    producing  a   i-hort- 
crested,    choppy  effect.      The   difficulty  arises   from  the 
fact   txiat   a   short-crested    sea   can   be    caused   by   one 
storm   system,    two    storm   systems    or   any  number   of    storm 
systems.      If  we    place    ourselves    on   a   snip   lying    to   in 
the    ocean,    maintaining  a    steady   true    comoass  heading, 
it    is    perfectly   possible    to   conceive    of   a   storm  system 
a  few  miles    off   the    bow   sending  wind  waves   toward    us 
and    another    storm  system  a   long  distance    off   the    star- 
board  quarter  sending   swell  waves    toward    us.      But    there 
could    just   as  well  have   been   just   one    storm   system  or 
more   than  two  and   they   could  have    been   off   in   many  com- 
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+>*0***H    TO    0**V€«T    €<*-*•    M    S*    S€A    $P*CT*A    TO    Lt>6    Sfc$*>fc    fOK»<    rf»N    »f 

1    FORMAT     ( 
. }    fOft»«AT    4i6f^4i 

3  FQkMAT     (7H1LAG    NO,6X<  a, 10*  i7A#llH    CORR    FT    M»4A 

4  FORMAT     <  2  *  .  13  »  3X  tS  (  l  1  <♦.«  ♦  < *  )  > 
4>lMfcN4U>*    (H)uU     ^AklJCl 

I • (^•m—w^X^9-  — •— 


f>0    22     1-1  tKAS^  S 

L     U     «     1.J0.10 
-      i2    -     I  1*9 

REA;>    2#  (CK(K  )  •    >k«  I  1  »I2  ) 

PRINT    3 

&©•   2?    J    *    *♦*© 

C    «     J«o 
AiC*C 
$I_«A«((°*C)***)»CMJ) 

QtViMI  <fl*€  >««frM.WkTE4UAiUJ-»  *• 

OM»B»C 

PHINT    4,     HCtC  .  ttF ( J  )  »5L tD£V 

2  3    CALL     SYSTEM 
6*C 


/&ufif  6 


n 


LAG  NO        OME^* 

o  .209UOOOOC  On 

7  .2U4"*oOOOP  0? 

b  .?79?0000r  00 

v  .3i4ioooor  00 

10  .TU9oO000r  On 

11  .383o0000p  On 

12  ,418poooop  oo 

13  .4537O00OF  00 

14  .48860000F  0* 
lb  .S23SO00OC  00 
lo  .558UO00OP  00 
17  .593««o000r  0*> 
lb  .ft28?O00OP  00 

19  •  <S631000'">F  On 
2u  .A98nO00OP  00 
21  .732QO000F  00 
2*  .767PO000F  On 

23  ,8027000^r  00 

24  .837*O00OF  On 
2b  •872FO00OP  On 

20  •907uO0OOe  00 
2  7  .Qtt2^0000P  On 
2b  .97790000^  00 
29  ,l01?l00nP  01 
3u  .1047000^  01 
31  ,l08l900*P  01 
3*  .lll(S8000P  01 
3 J  .H5l700nr  01 
3h  .118A6000P  01 
3;>  .1221500OF  01 


LN  OMFC-./* 
.1S635090F  OJ 
•1U0935A^F  01 
,1?7S82(S9F  01 
.115P0U3QF  01 
.ln52fc83UF  01 
•9S737318F  00 
.87036180F  00 
•7O031900F  00 
.7l62lli?F  00 
•6U7?182SF  00 

..5R2*797?F  00 
.5?205510F  00 

,.4*4^96*<3F  00 
.U10P29U7F  00 
.3*95361 &F  00 
.3107U601F  00 
•2*422600E  00 
#2i9774?3F  00 
•17721U*?F  00 
. 136392*2F  00 
.97171911E-01 
,5Q43l5fi^F-01 
.23063939F-01 
•12027381F-01 
,U«>9289-5^F-01 
.7ft7l87SSF-01 
.  H0«6745F  00 
.1U12391 IF  00 
•1710920PF  00 
.2000796 IF  00 


C0R9 
•55000000E 
•27000noOE 
.70000000E 
.60000000F 
•47000000E 
•  16700000F. 
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binations    of  distance    and   direction,    so    it    is    only- 
possible    to  express   what    is   happening  at   a   certain 
place    in   the    ocean   at   a   certain    time    by  a   forecasting 
or  hindcasting  technique.      Probably   the    most    typical 
situation   is   a   combination   of  wind   waves    caused   by  a 
local   storm  and    swell   coming  from  a  distant    storm  or 
storms.      This   is   a  matter  of   great   current   interest   to 
oceanographers   and    it   is   my  hope    that    in   another  five 
to   ten   years    they  will  have    accumulated   a   large   enough 
body  of  directional   sea   spectra   so  as   to  make   a  rea- 
sonably  sized    statistical  sample    possible. 

However,    for   the    moment,   we   must   content   our- 
selves  with  something   less   in   our   prediction  methods, 
namely,   either   the    use    of   long-crested    seas    or   "uni- 
storm"    short-crested    seas.      It   is   the    latter   that    I 
have   chosen.      Th"s   is   a   situation  where   one    storm  dis- 
turbance,   because    of   its   movement   and    irregular  winds 
causes   some    spread   in  the   wind  direction  and   thus    leads 
to  a  certain   type    of   short-cre stedness    in   the    seas. 
These    seas   are    characterized   by  a   large   amount   of  energy 
coming  from  a  certain  direction   and    progressively   lesser 
amounts   of  energy  as  we    look  away  from   that   direction 
until  there    is  essentially   zero  energy   coming  from  about 
70  degrees   on  either   side    of   the    original  direction. 
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This  directional  distribution   of   seaway 
energy  may   be    characterized    by  multiplying  a   long- 
crested   spectrum  by  a  "spreading  function*"     The   one 
I  have    chosen   is   based   on  the   work  of    Pier son,   Neumann 
and    James    In    (6)  and   is 

2  cos2//w 

rr 

The   result  of  this  directional  spreading  of 
seaway  energy  can  be   seen  at   the    top  of  Figure   9.      I 
have   chosen  this   technique   because    it    is   computationally 
practical  and    is   at   the   very   least,    one   step  closer  to 
reality  than  assuming    long-crested   seas*      It    is 
pertinent   to  note   that   my  experience  while  working 
at   the  David  Taylor  Model  Basin  showed   that   short- 
crested   predictions   by  this   method  are   consistently  8- 
9%  lower  than  long-crested  predictions   of  ship  motion; 
at    least   for  the   head   seas   case* 

With  the   foregoing  background   In  mind   it  is   in 
order  now  to  proceed  with  a  discussion  of   considerations 
that    led  to  a  final  choice   of   sea  spectra  to  represent 
the    sea   severity  range    of   interest   in   these    predictions* 

It  seems  reasonable   to  state   that  as  far  as 
ship  response    is   concerned,   there   are   two  characteris- 
tics  of  a   seaway  that  are    particularly   important,   namely 
the    total  energy,   which   is    proportional  to  the   area 
under  the    spectrum  representing   that    seaway,   and   the 
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distribution  of  energy  among   the   various   component 
waves,    an   indication   of  which   is    given   by   the    shape 
of   the    sea   spectrum.      One   neec    only   leaf   through    (7) 
to   see    that    nature    provides    us   with   sea   spectra   of 
many  different    sizes   and    shapes. 

Now,    if  we   wish   to    investigate    the   effect    of 
sea   severity   on   ship  response,    it    is    necessary   to 
establish  a    primary  weather   severity   parameter.      But 
in   doing  so,    caution  must    be   exercised    in    order   that 
other    parameters   affecting  ship-response   are    not 
ignored. 

1  have    chosen   significant   wave   height    (HW3) 
as    the    primary   sea   severity   parameter   because    it    is 
proportional   to   sea    spectrum  area  and    thus    total  sea- 
way energy.      (^1/3    is,    by  definition,    the   average 
height    of    the    one -third   highest   waves    in  a   seaway 
record.    (6)y  .      But,    in   order  not    to  neglect   the   effect 
of   sea    spectrum  shape    on    ship  response,    a  family   of 
sea   spectra  was   needed   where    at  each   increment    of   the 
primary   sea   severity   parameter,    Hwj,    there   would   be 
several   spectra   of   close    to   the    same    HWgj    but   having 
typical  vaiations    in   shape    among  them. 

The    I.Ioskowitz   family   of  mean   synoptic    sea 
spectra  found   in    (8)    is    the    best   family   presently 
available    most    closely   satisfying   the    above   requirements. 
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The  family  consists  of  five  mean  sea  spectra,  each 
spectrum  being  the  mean  of  8  to  14  component  spectra 
as  shown  in  Table  II.   See  Figure  8  for  a  plot  of  these 
five  mean  sea  spectra. 

TABLE  II 

Moskowitz  Synoptic*  Number   of  Component  Average   of  Component 
Mean  Spectrum  Name     Spectra  Spectra  H1/3   in.ft. 

20KT  12  7.4 

25  8  15 . 1 

30  12  20.4 

35  8  23.8 

40  14  32.2 

In  each  subset,    the   H3./3   of  the   component    spectra 
vary  in  most   cases   by  no  more    than   a   few   percent,    thus 
assuring  that   the    primary  differences   are   those   of 
spectral  shape. 

By  comparing  the   ship  response    caused   by   the 
Hwj   =  7.4   ft.   mean  spectrum  with   the   response   caused  by, 
say  the   H1/3   =  20.4   ft.   mean  spectrum,    one   can  determine 
the   effect   of   total  seaway  energy  on  response.      On  the 
other  hand,    by  comparing  the   ship  responses   due    to  each 
of   the    12   component   spectra   in   the   H-,  /*   ■  7.4   ft.    sub- 
set,   for  example,    one    is   able    to  evaluate    the  effect   of 
spectral  shape    on   ship  response    i.e.,   determine    the 


ft 
Moskowitz  classifies  his    spectra  on  the   basis   of  wind 

speed   in  knots. 
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variance    In  response  for  a  particular  Hwj,     To  my 
knowledge,    it  has  not  been  attempted   in  previous  work 
to  separate   these    two  effects   and    it    is   my  hope   that 
such  an  analysis   will  shed   some    light    in   our  quest  to 
learn  more  about    the  nature   of  ship  response   to  a  sea- 
way. 

The   Moskowitz  Synoptic  Family  of   Mean  Sea 
Spectra  converted  to  the    log-slope  form  is   shown  in 
Figure   8.      It    is  noted  that   Moskowitz  required  that 
the   sea  records   chosen  for  his  family  satisfy  certain 
synoptic   criteria  as  regards  non-transient  wind  condi- 
tions, duration,  fetch,  etc,  all  of  which  tend  to  insure 
that    the   seas  represented  are   fully  developed   or  close 
to  it.     He  has  also  eliminated   swell  from  all  records 
used.     So  one    must   keep   in  mind  that   the  ship  response 
predictions   correspond  to  the  nearly  fully  developed 
wind   sea  with  no  swell  situation  which,   of  course,    is 
not  a  precise  duplication  of  nature* 

When  there    is  available  a   larger  source   of   sea 
spectrum  data,    I  believe  the   soundest  approach  would 
be  to  pick  out  all  those   records  having  very  close   to  the 
same   H^/3  without   requiring  that   any  synoptic   criteria 
be   satisfied.     Thus,    if  the   sample  were    large  enough, 
one  would  be  reasonably  assured  that  a  representative 
variation  in  degree  of   sea  development  would  be 
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included.      Each  sample    of   records   having  the    same   Hwj, 
then,   would  have    the    variation   in   spectral   shape    pre- 
viously referred    to,    and    then   one    could    proceed   as    in 
this   analysis. 

Computation   of  Response  Data 

Once    the    prerequisite    information,   namely   sea 
spectra   and   RAO's    are    in  hand,    one   may   proceed   to   cal- 
culate   ship  response. 

A   verbal  description   of   the    computations   will 
be    given   to  accompany   the    graphic    illustration   given   in 
Figure    9.      First,    the    long-crested    sea    spectrum  is 
"spread"    or  directionalized   using   the   2   cos^/<w   spreading 
function   previously  described.      This    operation   results 
in  the    short-crested    sea   spectra   shown   by  dotted    lines. 

Next,    the   directional  RAO's   are    plotted    on   the 
same    Irw  abscissa   and   a  decision   is   made   as   to  wiiich  RAO's 
are   going   to   be   multiplied  by  which  sea  spectrum  component 
at  each   increment    of   In  oj  .      Figures    10  and    10a  will   aid 
in   this    "matching"    process.      Say,   for  example,    we   want 
to   compute    acceleration  response   for   the   head    seas   case 
i/1  -  180°)   as    is    shown  in  Figure    9.      We   align   the   heavy 
wind   arrow   on  Figure    10  with  the   Head   Seas   marker    on 
Figure    10a  and    see    that   the    180°  RAO  must   be    multiplied 
by  the    0°   sea    spectrum  component,    the    157   1/2/202    i/2° 
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RAO's    (which  are    the    same,    of   course,    by   symmetry) 
must    be    multiplied   by   the    22    1/2°   sea   spectrum  com- 
ponent   and    so   on   for   as   many   of   tne   RAO   components   as 
are    affected   by   the    various    sea   spectrum  components. 

This    is   done    at   each   increment    of    lnu  and   re- 
sults   in   the   nine   response    spectrum  components    (only 
five    curves    show   because    of   symmetry). 

Now   if   one    imagines   a   third   dimension   to   the 
Response    Spectrum  Comoonent    curves   whose    axis    is   normal 
to   the    plane    of   the    paper  and    whose    coordinate    is    /V, 
it    is    possible    to   visualize    the    integration  with  respect 
to  A  of   the    points   from   the   Response    Spectrum  Component 
curves   on   the    crossplot   for  each   In u>  increment   to  arrive 
at   a   point    on   the    Integrated  Response    Spectrum  for    that 
value    of   lnco  •      After   the    Integrated  Response   Spectrum 
is   thus   computed,    one   may   integrate    over  lnu  to   arrive 
at    the    area   under    it   which   is    the    value    of  ("u~E/ f  or    the 
case   considered,   namely  head    seas,    ship   length  =  370  ft. 
(\J"E   is   a   ship  response    parameter  that    is    discussed    in   de- 
tail  in  Appendix  A.      Briefly,    the  ~\jE   parameter  reflects 
certain  statistical   properties    of   a   short-term  sample 
record    of  response   during   invariant   conditions    of  weather, 
speed   and   heading.)     The   angular   integration   and    then   the 
In i>>  integration  were   cone    using  Simpson's   Rule    in   the 
computational   scheme   employed.      Nineteen   In  u)  increments 
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of   0.1  each  were    used   along  with   the    22   l/2°  angular 
increments   as    shown.      The    right    and    left  hand   extremities 
of    the    lnv*3scale,    namely,    -1.6   and    +0.2,    respectively, 
were   determined   by  the    termination   of  most    typical  sea 
spectra   on   the   right  hand   end   and    the    termination   of    the 
RAO's    on   the    left  hand   end. 

The    process    described    above    can   then   be    done    for 
the   eight    other  headings    between  head    and    following   seas 
so   that  each   sea   spectrum  and   each   set   of  RAO's   for  a 
given   ship   len    th  results    in  nine  IfE  values,    one    for 
each  heading  from  head    to   following   seas   at  22    l/2° 
increments . 

This   was   done   for   the    five    Moskowitz  Synoptic 
Mean  Spectra  for   each   of   the    five    ship   lengths   and   at 
nine   headings,    resulting   in   the   225  "\jE  values   found    in 
Figure    11. 

Since    each   one    of    these  ME   values    is    the    result 
of  the    process    of   commutation   just   described,  where    sea 
spectrum  components   are   multiplied   by  RAO   components    to 
arrive    at   the   Response    Spectrum  Components   which  are    in- 
tegrated   to   arrive    at    the    Integrated   Response    Spectrum 
which  in    turn    is    integrated    to  arrive    at    the  \fE   value, 
it   can  readily   be    seen   that   a   very   large   number   of 
numerical  calculations    are    involved.      However,    since 
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PIGIHE   11 
!£  MEAN   1^  ACCELERATION    /ALUES    (x)   in  g's 

^l/5   =  7»4'( Moakowitz  20KT  Synoptic   Mean  Spectrum) 


Ship 

Length    — 

—>     248' 

3  02' 

370' 

452' 

552' 

180° 

.2108 

.1676 

.1238 

.0885 

.0594 

157   1/2 

.2021 

.1631 

.1212 

.0882 

.0594 

135 

.1733 

.1476 

.1148 

.0374 

.0625 

112    1/2 

.1467 

.1261 

.0999 

.0776 

.0559 

90 

.1037 

.0946 

.0812 

.0676 

.0516 

67   1/2 

.06  GO 

.0621 

.0544 

.0453 

.0348 

45 

.0352 

.0354 

.0335 

.0296 

.0233 

22    1/2 

.0157 

.0158 

.0151 

.0130 

.0101 

0 

.0090 

.0090 

.0090 

.0085 

.0063 
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W3   =  13. l'  (Moskowitz  25KT  Synoptic   Mean  Spectrum) 


130° 

.3181 

.2843 

.2467 

.2062 

.1623 

157    1/2 

.3033 

.2727 

.2369 

.1993 

.1574 

135 

.2644 

.2408 

.2122 

.1815 

.14  62 

112    1/2 

.2117 

.1967 

.1747 

.1521 

.1242 

90 

.1470 

.1412 

.1298 

.1168 

.0985 

67    1/2 

.0921 

.0901 

.0338 

.0758 

.0646 

45 

.0489 

.0497 

.0484 

.0452 

.0393 

22    1/2 

.0222 

.0224 

.0218 

.0199 

.0173 

0 

.0123 

.0128 

.0128 

.0124 

.0105 

Hl/3   " 

20. 4 '    (Mo3kow 

itz  30KT 

Synoptic 

Mean  Spec 

trura) 

180° 

.3753 

.3447 

.3106 

.2742 

.2356 

157    1/2 

.3576 

,3300 

.2975 

.2634 

.2264 

135 

.3112 

.2901 

.2641 

.2363 

.2047 

112    1/2 

.2481 

.2349 

.2147 

.1942 

.1689 

90 

.1722 

.1675 

.1572 

.14:3 

.1283 

67    1/2 

.1079 

.1066 

.1008 

.0934 

.0828 

45 

.0578 

.0588 

.0576 

.0547 

.0489 

22   1/2 

.0265 

.0267 

.0262 

.0244 

.0217 

0 

.0154 

.0154 

.0154 

.0150 

.0131 
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FIGURE   11    (Cont'd) 


SAMPLE  MEAN    TfB  ACCELERATION    TALUKS    (5t)    in  g's 
Hw3   =  25.8       (:.iOskowitz  35KT  Synoptic    Mean   Spectrum) 


Ship 
Length 

{ 

180° 

157    1/2 

135 

112    1/2 

90 

67    1/2 

45 

22    1/2 

0 


248« 


302' 


370' 


452' 


552' 


.4287 

.3994 

.3646 

.3252 

.2823 

,4082 

,3819 

,5488 

.3121 

,2711 

.3548 

.3348 

.3085 

.2787 

.2439 

.2817 

.2698 

.2496 

.2282 

.2006 

.1951 

.1914 

.1812 

.1691 

.1509 

.1222 

.1214 

.  11<^8 

.1084 

.0972 

.0655 

.0667 

.0557 

.0628 

.0569 

.0302 

.0304 

.0299 

.0280 

.0254 

.0176 

.0176 

.0176 

.0171 

.0152 

^l/3   =  32. 21      (I'.loskowitzi  40  KT  Synoptic   Mean  Spectrum) 


180 

157    1/2 

135 

112    1/2 

90 

67    1/2 

45 

22    1/2 

0 


.4732 

.4524 

,4220 

.3871 

.3485 

.4554 

.4323 

.4032 

.3707 

.3337 

.3957 

.3734 

.3554 

.3239 

.2973 

.3139 

.3038 

.2358 

.2668 

.2416 

.2179 

.2153 

.2064 

.1955 

.1786 

.1369 

.1366 

.  13 15 

.1248 

.1142 

.0741 

.0754 

.0745 

.0713 

.0660 

.0345 

.0347 

.0341 

.0323 

.0296 

.0201 

.0201 

.0201 

.0196 

.0177 
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these   calculations  are   of  a  repetitive  nature ,    including 
several  cyclic   permutations,    they  are   especially  amenable 
to  calculation  by  electronic   computer,   and   this    is   how 
they  were    carried   out. 

The   program  shown  in  Figure    12,  based   on  a 
similar  one    in  use  at   Webb   Institute,  was  written  by  the 
author  in  the  Fortran   II  source    language  while   on  temporary 
duty  at   the  David  Taylor  Model  Basin  and  the   calculations 
were   carried   out   on  the   UNIVAC   IARC   facility  there •     A 
sample   of  the   print-out   can    be   seen  in  Figure    13,      It 
might   interest   the   reader  to  know  that  the   computer  took 
less   than  three  minutes   to   compute   the   225  TJE  values   and 
having  done   similar  calculations  with  a  desk  computer, 
I  estimate    it  would   take   about  a   month  of   8-hour  days 
to  calculate   the    same  number   of  iS  values   manually. 
Thus    it    is   seen  how  the  electronic   computer  can    greatly 
enhance   the   scope   of  an   investigation  like   this   one. 

The  tabulation  of  ~1JB  values  in  Figure  11  shows 
the  effect  of  seaway  total  energy  on  acceleration  res- 
ponse, but  it  does  not  show  the  effect  of  sea  spectrum 
shape* 

As   previously  described,    the  effect   of   sea  spectrum 
shape   on  ship  response  was   to  be   evaluated   by  computing 
"\fE  values  for  each  of   the   component   sea  spectra  in  the 
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LABEL       CKSRD 

C      R5C  IS  RESPONSE  SPECTRUM  COMPONENT.  S  IS  SEA  SPECTNU* 
-« RSI  !3  iNtCORATEt)  l>Cy^XV£(^OM,  5*-+$--^^****^-  >»j^f  j pWf<* 

C      RAO  IS  RESPONSE  OPERATOR.  S^t        Is  AGAINST  LN  0«E<»A 

*R€A  t*,r>€-R  «-S+  fc*t*M.-$  0.e»*3H-SO*>,  4  IS  R*4  **  64*&4.-£  A*»4l 
1  FORMAT  (6F7.4) 
?  FO*a**T-  fi-GE?.** 
3  FORMAT  Mi) 

^-rOft*UT     f  M6»  SM1P    LENGTH     -f  

P     FORv*T     OOM'Rvs    KfSPCNbE     FOR    *U     SU8    0    UF     .F6.1.4H     IS     .£12, 
V    FORMAT     -fh*l,    47>f*      13H     S*A     SPfc<IT«<*4> 

DIMENSION     SHII9I.     RSK1V).     X  I  1  7  .  1  9  .  5  >  .    KCI17),     PR  (  1  7  )  .     S(19l        Tf 

-A    -     fr^-^?4»  I 

3 


8    ■     C.262 


5M     (?)  ■     <♦  . 

**    *3f  -2. 

*u     (<♦)  ■     A. 

S*4    <5l  •    2, 

SM     (6)  «A. 

*M     { 8  )     •     4. 
SM( 1C)     «     A. 

***t*H    •-*• 

SM(  ID      -4. 
*«<  1»>    — -g-»- 


■•"* 


VMU)  «  <♦. 
S'M  16  >  *  4. 
SM( 18)     •    4. 


1 

3 


j 

1 


&O    3  9    L    *     1.5 

READ    l.XIl  •I»L>»Xi2tItL)»X(3,ltL>tAU»ItLI»XC5tI»L>tXI4»ltU 

1»X(13,I,L>.,.  (l4.I,L).*(15»I.L».X(l6.ItL).A(17»I.L)  -;  < 


KEAO 


03»    CQ*TlN«>g 

READ    J.    KAStb 
ttO- 4^    f    •     }r*A&£* 


1 


^A*   ItWUUMIJWHl^tSiUI.Hl^^^^^l?  f  »4444*  ♦**  *** 
PR  INT    9 

♦O    49    L    «>-4-r*- 


PRINT     7  "^ 

OMtrtr  •    d»^*«f7  

DO    6  9    H\JZ     .     1,9 

S4#*     *     £w& 

PMUZ     »    MUZ- 1  """* 

DO    64    10    •      !.!<?ti  — — f 
I-4i«f    •    £<•© 


LLO    ■    LC-    3  ~> 

AMUW     *     -1.5708  ~* 

RC(JMU)     .     X(  J^'J,L0»L)»S(lC)»0.637MCoVf(AHUW)  )»#2  """"^ 
A«(^44>M*R-C  iUM^»«4 

Fo*Te»NJ[  Ceurcr  fo***„  t0  Colore  Acrt**Ar,*i  &s+~,t*r 

3£  F****/2 


\ 


PJMU-JMO-I 


RSIIL0I-R5H  JMU! 

JMK* JMU*1 
Ofc*    £0*Tl*ttf 

SUM«l75l  (LO)»SM(LU)*bUM 

O«5QRTF(SUH*C03  3n) 

-P*4**  *^PW**^ 

069    CONTINUE 

END 
Mjlll 


^a 


i 
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F/&v*£  /2  Ct^TP 


sitt  «53£l;tR(ji 


ftCUtLtK*TIO^  Ktar'O^st 


SHIP    LCN*TM          1                ■*f6 

J.J 
2?.* 

4  5.0 

57.* 

?0.0 

IMO 
157.5 
1  95.0 

lb 
TS 

TS 

?s 

IS 
TS 
'S 
TS 
IS 

*M5 

•  MS 
*M5 

•  MS 

•  MS 

•  MS 

PR 

•  MS 

•  MS 

^E«PUN>*    po<?    MJ    bJ*    0   OF" 
JPSPONS*    FM    MJ    SJ5    0    OF 

•  re;poNSr    F**    MJ    SJ*    0    of 
^rsPON«»f   P->9.    MJ    SJ«    0    OF 
•rsPONS*    F*9    MJ    SJ*    0    of 
*ESPON*c    F*e*    Ml    SJ*    0    MP 

•  <rs»0rt»*    P**    MJ    SJ»    0    oe 
^rsPONJsr    f*9    MJ    $j*    n    op 
^rsP0N"sP    P*P,    MJ    SJa    0    oP 

J.?3l30£-J2 
3.1S3-»E-31 
0.3?0-«P-01 

0.«,33iE-J1 
0.109-»E    JO 
J.1S21E    JO 
3.1^2'E    JO 
3.>0SoE    00 
0.?l'4iE    00 

•  MS 

SHIP    LtNuiTu         i 

^rSPONsf    P«r*    MJ     J>J-«     0    QP 

3. J 

IS 

J.?3SO£-J^ 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•.CAPONS*     F»H     MJ     vj*     o    OP 
•PSPON^e    F-P.    MJ    SJ*    0    of 
•esPONsr    p«n    MJ    Sj*    o    of 

•  'SPON-sr    F*9    MJ    SJ*    0    OF 
•FSPONS?    ?«•    MJ    SJ'    0    OF 

•  eSPON-se    F«n    MJ    SJ*    f)    OF 

4  3.  0 

sr.s 
»oo 

1  1  ?.«s 
1  *3.0 

TS 
TS 
TS 
TS 
TS 

0.1S"4«.E-J1 
3.**>«iiE-01 
3.5371E-01 
0.9*,2-»E-01 
0.12SCE    JO 
J.l'46'iE    00 

•  MS 

•  MS 

9-S°0N«»r    F*P.    MJ    tj«j    0    *>c 
•^S^ONS*    F-p    MJ    ?.J*    f)    of 

SHIP    LtN»m          *     - 
^ESPON«»ff    F«n«    MJ    SJ-9    0    OP 
•^SPONsr    p.p.    MJ    sj*    0    ** 

MN1 
i  93.0 

TS 
TS 

IS 
TS 

1.160'tE    JO 
J.1S4-»E    30 

•  MS 

•  MS 

>.J 
2?.S 

J.9330E-J2 
0.157-E-Jl 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

^<BONSf     F*H     M  1     SJ»     0    *c 
^=-SPON*r    F*p    MJ    SJ*    0    OF 
VCUNsr    F">3    MJ    SJ*    0    or  . 
•.-SPONsr    F«p    MJ    *,ja    0    nF 
•.ESPONsff    F*p    MJ    SJ*    O    oP 
^cSPONse    F-*.     Ml     t,J"     0    OP 

45.3 

sr.s 

135.0 
157.5 

TS 

TS 

TS 
TS 

TS 

tS 

3.3'49-»E-01 
0.3<59"«E-31 
0.333*>E-0I 
0.100«E    JO 
0.11<4<>E    J1 
0.12U^E    00 

•  MS 

•.FSPON'sF    p*p    MJ    SJ»    0    OP" 

19}. 0 

TS 

0.121->E    00 

*M$ 

•  MS 

•  MS 

SHIP    LtNaTu          4 
pr«.pgNs»    F«R    MJ    *,J9    0    OP 
•csPONiF    ?«•    MJ    SJ*    0    OP 
•  FSPONs«    F«p    MJ    Sja    o    of 

3. J 
2>.5 

45.0 

IS 
TS 

TS 

J.97>4^E-J2 
0.134'E-01 
J.105OE-01 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  ^SBCJNbc    F»P    MU     S>J*    0    «F 
•.fspjNs*    p»*    Ml    sj=<    o    OF 
•rsPONi?    P«R    MJ    sjs    o    OP 
•.FSPONarr    P-P.    MJ    SJ*    0    OF 

•  PAPONS*    F«R    *J    SJ<9    TO    OF 
^-SPONSr    F*9    MJ    Sja    o    of 

ST.s 

J-J.O 

112. S 

195.0 
157.5 
1  9  ■)  .  > 

TS 
TS 

TS 
TS 

TS 
TS 

D.'4^J->E-01 
O.S92->E-01 
0.»63^E-01 
0.^52»?-01 
O.S35«E-0! 
J."'96-'E-01 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 

•  MS 


SHIP  UtSaTu  ■> 

•  ESPUNPF  F->S»  MJ 
•,F«tPJN!»»  P-i9  MJ 
RESPONSE  F->H  MJ 
•.ESOONsc  F«p  MJ 
•.ESPONsff  F^H  MJ 

•  ESPONsf  F->R  MJ 
^ESPONsr  F«,p  MJ 
•ESPONpf  F»R  MJ 
9E^»0N!»*  P«9  MJ 


bJ9  O  Of 

Sja,  O  of 

SJ*  O  or 

SJ*  0  op 


SJ*  o  or 

SJ^  0  OF 

SJ»  0  OF 

SJ'  0  OP 


O.J  IS 

22.5  TS 

4  3.0  TS 

47.5  TS 


112.5 
135.0 
157.5 
13->.0 


TS 

TS 

TS 

TS 


J.SS«i-»E-J2 
O.IOPE-Jl 
J.?^'4Ti;-Jl 
3.3'41-E-Ol 
D.50S-»E-0l 
0.522*E-01 
O.STZ'JE-Ol 

0.553'E-Ol 
J.4?a«E.^j 
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five    Moskowitz    subsets.      Since    there    are    a  total   of   54 
component    spectra    in    the    Moskowitz   Synoptic    Family   of   five 
mean    ?pectra,    converting  each   of    these    to    the    log-slope 
form,    plotting   them,    tabulating   and    card-punching  them 
as    input   to   the    ship  response    computer   program  would 
have    been    too  formidable    a   task  to   undertake    in   this 
thesis.      Instead,    the   above   was   done    only   for   the    12 
component    spectra   of   the    H3/3    =  7.4    ft.    mean   spectrum 
and   the    14    cornoonent    spectra   of   the    HW3   =  32.2  ft.    mean 
spectrum   so   as    to  bracket    the    range    of   the    primary   sea 
severity   parameter  H3/3 .      The    response    characteristics 
for   intermediate    values    of  H3/3   were    obtained   using 
straight    line    interpolation   which   turned    out    to  be    a 
quite   reasonable    procedure    as   will   be    discussed    later. 

The  *\M  values   for   all  nine   headings   and    five 
ship   lengths    in  the    12   component    sea   spectra    of   the 
Ht/3   =  7.4   ft.    subset   and    in   the    14   component    sea   spectra 
of   the    Hn  a?    =  32.2   ft.    subset   were    computed    using    the 
same    computational   scheme    as   was    previously  described    and 
then   for  each   combination   of    ship   length  and  heading,    the 
following    important   characteristics    of   the  TJE  distribu- 
tion were    determined: 

(a)      the    sample   mean   of  "\fE  was   computed   using 


/n 


X     =     1        <C       Xi 

n     £r, 


where   n  =  12   or   14  depending   on  which  subset 
was   under   consideration  and  x  =  sample   mean 
1/1  and   xi   *    ("tflJi 
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(b)   the  variance  of  l/E  was  computed  using: 
s2=  1  (I   Ui)2  -   n  A 


again  where  n  =  12  or  14  and  s2  =  sample 
variance  and  s  =  sample  deviation 

These  computations  also  were  done  on  the  DTMB 
UNIVAC  LARC  using  the  Fortran  II  source  program  written 
for  this  purpose  which  is  shown  in  Figure  14.   This 
program  is  similar  to  the  one  in  Figure  12  except  that 
a  section  for  computation  of  sample  means  and  devia- 
tions was  added  onto  the  section  that  computed  the  uE 
values.   The  print -out  of  results  from  this  program 
may  be  seen  in  Figures  15  and  15a. 

At  this  point,  one  small  issue  was  disposed  of. 
For  each  ship  length  and  heading,  the  mean  'tfE  values  in 

Figure  15,  for  example,  should  agree  with  the  llE  values 

»/■■ 
given  in  Figure  11  for  the  H3/3  =  7.4  ft.  mean  sea 

spectrum,  that  is  to  say  that  the  mean  value  of  ^E 

obtained  using  the  subset  component  spectra  should  be 

equal  to  the  1JE  obtained  using  the  mean  sea  spectrum. 

This  is  mainly  a  crosscheck  on  the  computational  scheme 

used.   A  comparison  of  the  data  shows  a  maximum  difference 

in  the  order  of  4  percent,  so  the  procedure  of  using  a 

mean  sea  spectrum  to  obtain  mean  IfE  response  seems 

justified.   This  conclusion  engenders  confidence  in  the 

*\$E  values  shown  in  Figure  11,  because  the  ones  for  the 
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FOh  UP'  HI  vTA    PL  >■  ■  r\_  i 

bo 

•  • 

:     LN 
SPECT  ft    S    IS  PECTRUM 

-  E  A  :-...,  .         •  »  . 

^fi     Is    SAMPLt     MEA.N 

1  FORMAT     f6F7*41 

2  J  .2  .' 

i  ■ 
7    •  '     (  /]  SHIP    tt.V    T  -.     ♦  I  3  ) 

»F6«1»4H  »E1  2.  -4  i 

»    37X» ' 13H  •  ..  •  X »22H    ACCELERAT I 

4F     Unli  •  .  <»  '  18ri     VLLUC! 

■     hi  •  -       <* . '  rn     •  20x»22i      i     -Actwi  n    *es 

>     42X»     24H  •  IS) 

31  »    44X*    2  OH    VELOCITY    C 

36  AT    (34  EVIAT3  Q    Oi      »F6«1«4H     IS    »Ei4.8»13 

i     14.8) 
38  '  I  »     ■*'•  X  ♦     24H    D  i  /J  AT  IONS) 

■  •.      >  :  »  v  j  ,     <  (  1  7  ,  1 9  *  5  )  i     3 Q  ( 1 7  J  1    PR ( 1 7 ) »     5  f  1 9  )  ♦    Q  (  26 

... 

•    -    S.262 

■■    ( l )    =    :  . 

( 3 )  -     . 

SW     (  4  )  *    <•  . 

SM    15)  =    ?. 

!  6  )  =  *. 

'     {  7 )  s     /  . 

•     (  e  )  =    4  • 

5M     (Qi  *  2. 

SM  (  1  0  )  s  4  « 

5M( 1  1  )  =  ?, 

■  (  I  ?  )  -  •+  . 

1(13)  «  • 

M 14)     =    4. 

SM  (  I  5  )     =    2 . 
1(16)     = 

t :  7  i    s   2  . 

I U  8  J     s    ^  . 
SMU9.I      a     i. 
DC    1 9    L.l    =    1  *  3 
DO  a     U5 

-J  9       I       r      1,1. 

READ    1  »X  ( I » I  •  L  )  »X  (  2  1 1  »L  )  *  X  (  3  *  I  f.L  )  •  X  {  4»  I  » L )  »X  (  5  »  It  L)  »X '(  6  »  I  >L  J 

,        READ     1  ,  X  <  7  ,  I  ,  L  )  »  X  .'  d  ,  I  .  L  )  » X  {  9  *  I  *  L  )  *  X  i  i  0  ,  I  ,  L  )  »  X  (  i  1  ♦  i  »  L  )  »  X  (  1  2  ♦  1  »  L  ) 
READ     1  *  X  (  i  3  *  I  .  L)  »  X  (  1  4  ,  I  »  l  )  »  X  (  i  b  »  I  ♦  L  i  *  X  (  I .6  »  I  »  L)  »  X  (  I  ,'  ,  I  »  L  ) 
039    CONTINUE 

READ    3»    <AS- 
DO    6  9     I     =     1»KASES 
3  2    READ    2  *  S  (  1  >  ♦  S  (  k  )  ♦  5  {  3  >  » S  (  4  )  *  b  <  :>  )  *  S  (  6  )  *  S  I  7. }  » S  C  8  )  » S  (  9  j  »  3  <  10  ) 

AD    2*S<  11)  »S(  12)  »3U3J  »S(  14)  jS(15)  »S<i6) »S(  17)  »6<  18)  »S(  19  J 
(33»34»35) •    LI 
P  R  I  N  T     9 
TD    37 
►    PRENT    ir 
GO    TO    37 

For  TA^'IT  jQuecjF -  Pp£6f?AM   TO  CcfrpVTJE-    ^AM/V..rr  fUr/*JS  }  T^JE-WAT'dMS 


r  n 

*  I    ■;■   i    69    L     -     ll 
I     7,    L 

jZ-s    1»9 

-  1 

JZ 

t «    ?  t  i 

-  3 
j  M  U     - 

=    -  1 . 

»5»] 

•  ...  ■.  )  )  **2 

-  i    I  L       )  H 

■ 

'  ♦  5  7 1      V      •     } )  •  N 
57    RC(J  -      ■  ■  »LO»LH  »63  7  -  2 

PJMU     =     ..'■'    -  i 

.--:.    *   p  j  ■•      -   9  c .  o 

[  (LOJ=i 

-  ! 
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<    < 

XK     = 

r  C     (  b  5  •  5  6  %     3  j  >    L  i 

:  R  I  \i  T     30 

ro  \2 

56    PRINT 

:    TO  12 

PRINT  J8 

12    00     16  L     ■     i  » :> 

PRINT      7,      |. 

SUM     =     0.0 

rsuM  =   • . 

1?    DO     15     I     =     1  ,  ■ 

I  I  tL  >y,JZ  )     '-    Q(  I  ,l  »MUZ  >**2 

SUM    -     SUM    +    Q( I »L -MUZ ) 

TSUM  a  TSU,M  +  QS(  I  »L»KU2  ) 
15  CONTINUE 

XB  =  SUM/XK 

R  =  XK* ( { XB)**2 ) 

VAR  =  (  1  ,0/  (  XK-:  .C •)  )*  (  TSUW  -  RJ 

DEV  =  SQRTF(VAR) 

Nl     =    MUZ    -    1 

Nl    a    22f5*f|] 

PRINT    36 »    Nit     DEV,     XP 
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16    CONTIN 

*    < 
X  K     =     ' 
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GO    TO    8  8 
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r 
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L  7    I    a 
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three   intermediate   values   of  H^/g  were   computed   using   the 
mean  sea   spectra   only. 

The    sample  i|E  deviations  were    linearly  inter- 
polated  over   the   Hw^  range   and   are   shown   in  Figures 
16  a-e   for   the   various   ship   lengths   and  headings. 

To  summarize,    it    is   seen  that     at   this    point, 
enough   information   is   available   so  that   the  effect    on 
acceleration  response   of  both  sea  spectrum  total  energy 
and   sea   spectrum  shape   can  be   examined. 

The  effect   of  »a  variation   in  sea   spectrum  total 
energy  manifests   itself    in   the   five   groups   of   sample 
mean  "\$E  values   in  Figure    11.     The   effect   of  a  variation 
in   sea  spectrum  shape   manifests   itself   in   the   sample 

uB  deviations   shown   in  Figures   16  a-e;   these   deviations 
from  the   mean  exist   solely  because   of  this   shape  effect 
on   ship  response. 

At    this   point,   the  ill  values   in  Figure   11  could 
be  variously  plotted  to  show  the   trends   of  1$E  stern 
acceleration  with  heading,    length  and    sea   severity,   but 
before   this    is   done,    it    is  desirable    to  treat   these   IJB 
values  which  are   characteristic   of   short-term  response, 
as  a   statistical  sample   and   use   them  to  achieve   a  rea- 
sonable  estimate   of  what   $B  values   might   be  expected  within 
the   total  universe   or  population  of  all  possible  destroyer 
stern  acceleration  ~vJE  values  during  a  destroyer's    lifetime 
under  similar   conditions   of  heading,    length  and   sea   severity. 
The  next  section  will  be  devoted   to  this  estimation. 
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STATISTICAL  ESTIMATION   OF    POPULATION  ^  VALUES 

As  discussed    in    the    previous    section,    at  each 
combination   of  heaaing,ship   length   and   Hwg,    there    is    a 
mean  value    of  ]JE   (Figure    11)   and   a  distribution   of  1JE 
values   above   and    belov/  that   mean,   which  distribution    is 
characterized   by  the    sample   deviation  found   in  Figures 
16  a-e    for    the    same    combination   of  heading,    length  and 

Hl/3» 

The  assumption  was  made  that  this  distribution 

is  Gaussian  or  normal.   The  normal  probability  distribu- 
tion is  completely  defined  by  a  mean  value  and  the  de- 
viation from  that  mean.   The  primary  defense  of  this 
assumption  of  normality  is  that  it  does  not  seem  un- 
reasonable in  the  light  of  present  knowledge*   Perhaps 
when  much  more  work  in  this  area  has  been  done,  a  slightly 
better  distribution  can  be  found.   The  size  of  the  statis- 
tical sample  used  here  is  too  small  to  permit  construction 
of  a  valid  histogram.   This  assumption  of  normality  has 
also  been  used  by  other  investigators  (12). 

The  next  and  last  hurdle  is  how  to  use  a  sample 
mean  and  deviation  to  predict  the  true  total  population 
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mean  'tfE  and    its   deviation.      This    is   where    one   may  make 
good   use    of   the    concept   of  statistical   inference,    i.e., 
one   may   infer  from  the    characteristics   of  a   statistical 
sample    the    characteristics    of    the    total  population.      In 
this   case    the    assumption  was   made    that    in   the    sample, 
the  \fE  values   are   normally  distributed  with  a  mean,? 
and   a  deviation,    s. 

Sample :      N    (x,a2  ) 
In  the   population,    the  "|/E  values  are   also  assumed 
normally  distributed,   with  a  meany*<and   a  deviation,    (T. 
Populat  ion :     N ( /f , <T2  ) 

The    task  now  is   to  estimate/?  and  <3~ on  the 
basis   of  2  and   s.      Presently  available    theoretical 
statistics   gives   one    two  choices   that   are   quite   differ- 
ent     in  concept,   namely  the    processes   of  point  esti- 
mation and   interval  estimation. 

I  will  first  discuss   point  estimation.      Maximum 
likelihood  estimators  may  be   derived  which  give   a  con- 
sistent,unbiased  estimate   of  the    population   parameters. 
In  the   case   under  consideration,    the    unbiased  maximum 
likelihood  estimator  forAis: 

A»  i   £**    =  x  (9) 

n 

and   the   unbiased  maximum  likelihood  estimator  for  tf~*is: 


<f  =  )|  <    [Xk  -X)8      -  s 


&i1*' 


*It   can  be    shown   that   this  expression   is  equivalent   to 
the    one   used   on   page 40. 
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This  means  that  the  best  estimate,  one  can  make 
in  this  case,  of  /4  and  (f  is  to  use  £"  and  s,  respectively. 
However  this  somewhat  vague  approach  is  not  very  attrac- 
tive unless  the  sample  is  of  a  respectable  size,  which 
in  this  investigation  it  is  not.   Perhaps  the  largest 

drawback  to  the  use  of  point  estimation  is  that  when 

I 
sample  size  is  relatively  small,  one  is  provided  with 

no  specific  expression  of  confidence  that  can  be  had  in 

the  estimation  of  the  true  parameters  H  and  (f • 

It  is  due  to  the  above  that  I  decided  to  use  the 
concept  of  interval  estimation.  With  a  specific  level 
of  confidence  in  mind,  say  90%,  one  may  predict  a  para- 
meter interval  or  region  which,  withM  90%  confidence, 
will  cover  the  true  parameter  value  or  values  being 
s  ought • 

It  has  been  shown  (9)  that  the  statistic 

<ro/iTn 

is  normally  distributed  with  zero  mean  and  unit  variance, 

and  that  the  statistic 

ns2  - 

752-  * 

is  chi-square  distributed,   where   (To  and  A>  are   the    true 
population  parameter  values. 

Thus,   using  any  convenient   set   of   probability 
tables,    the    statement  may  be   made   that 

*    These  two  statistics  are  stochastically  independent. 
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Pr   F-1.15   /    x  ->o    /    1.15        ■   .875 

t  7o7vr  *■       J 

from  the   standardized  normal  distribution  and   that 
Pr       4.57    /      ns2      /    19.7    |        -   .900 

from  the  y^  distribution  with  n-1  degrees   of  freedom 

where  n  =  12   in  this  case,    (n  =»  12  was   used  for  the   sake 

of  uniformity  even  though  the  H1/3   ■  32.2   subset  had   14 

component  spectra-this  will  yield   slightly  conservative 

results ) 

The    joint  probability  is  (9) 

P   f-1.15^  X  ~/r0    /   1.15,   4.57  ^  ns2    ^19.7     = 
[  CTo/Yn'     >  0*0%  J 

(.875) (900)    =    .788 

These   inequalities  define  a  confidence  region   in  the 

tf  9/Lf  parameter   space   as   shown   in  Figure   17. 

There    is   a  78,8%  chance   that   the   crosshatched 
parameter  region  covers   the   true    parameter   point    (A>,$o), 
wherever   it  may  be,    that   is   to   say   that   if  100  samples   of 
12  each  were   chosen  and   if  the    100  resulting  78.8%  con- 
fidence  regions  were   plotted  as   above,   about   79   of  them 
would  cover   the    true    parameter   point   and   the   remaining 
21  would  fail  to  do  so.     This   is   the   quintessence   of 
the   concept   of   interval  estimation. 

The   approximately  80%  confidence    level  was   chosen 
as   being  reasonable   and  the   following  discussion   is 
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intended  to  show  that  it  is  more  than  conservative 
enough  for  the  purpose  of  these  estimations, 

The  maximum  likelihood  estimators  tell  us  that 
P  is  the  most  likely  location  of  (A>,  (To),  With  this 
and  the  approximately  80%  confidence  level  in  mind,  I 
have  chosen  (Mi,  6  u)  which  will  give  the  "worst11  normal 
distribution  in  terms  of  extreme  values  as  the  true 
parameter  point.   This  is  a  fairly  conservative  choice, 
beoause  of  the  likelihood  that  (  <f  o,A?)  is  somewhere 
near  P.  I  feel  this  is  a  situation  where  a  little 
judicious  conservatism  is  better  than  being  inordinately 
bold  in  ignorance • 

The  parameters  of  N(Ai,  Cvr)   may  thus  be  evalua- 
ted from  the  following  simple  relationships: 
Ai  =  x  +  .540  s 
fa   =  1.620  s 

These  equations  show  the  amount  by  which  the 
sample  parameters  are  increased  so  as  to  arrive  at  an 
estimate  of  the  population  parameters. 
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INVESTIGATION   OP   TRENDS   IN  iJf  ACCELERATION 

Prom   this    point,    mean  l|E  acceleration  will   imply 
estimated   population     mean  1JE  acceleration    (A*u)   and   vE 
acceleration  deviation  will  imply  estimated    population 
deviation   of  "\JE   acceleration     (  0  u). 

From  the   data   in  Figure   11    (x)   and  Figures   16 
a-e    (s),    many  combinations   of   ship  heading,    length  and 
**l/3   can  be    vised   to  investigate    trends   in  *\f§  accelera- 
tion.     For   purposes   of   the    present  work,    I  have   considered: 

(a)      one    ship  length   (370  ft)   -  all  headings- 
all  S1/3 

(b)      one   heading   {head   seas)    -  all  lengths- 
all  Hl/5 

These   cases   show   the   effect   of   ship  length  and  heading 

as  well  as  H-^  /^   on  mean  ~\B  acceleration  response  and   its 

deviation.      Other  cases   could,    of  course,    be   considered 

if  desired,   using  other  data  from  Figures   11  and   16  a-e. 

Figures   18  and   18a  deal  with  the  370-foot   ship 

length  case   and   show  how  mean  "\JE  acceleration  and   its 

deviation  vary  with  heading  and  Hw^.      Figures   19  and   19a 

show  trends  with   ship   length  and   H3/3   for   the   head   seas 

case. 
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Finally,  by  crossplotting  Figures  18  and  18a 
at  discrete  values  of  H1/3  (which  can  be  thought  of  as 
different  sea  states  -  see  Appendix  B),  Figure  20  was 
derived.   Similarly  crossplotting  Figures  19  and  19a  led 
to  Figure  21.   In  the  interest  of  clarity,  only  enough 
deviations  are  plotted  on  Figure  20  to  demonstrate  the 
fact  that  the  deviations  hehave  in  much  the  same  fashion 
as  do  the  mean  'M   values  for  any  given  ship  length,  i.e., 
a  plot  of  deviations  results  in  curves  similar  to  the 
curves  of  mean  TJE  acceleration  in  Figure  20. 

If  one  accepts  this  conclusion  of  similar  be- 
havior of  means  and  deviations,  the  fairly  mild  curva- 
ture of  the  mean  "\JE  curves  between  H3/3  =  7.4  and  32,2 
feet  in  Figure  18  tends  to  justify  the  assumption  of 
linearity  used  in  interpolating  the  sample  deviations  in 
Figures  16  a-e  in  the  same  Hw,  range. 

There  are  some  important  points  to  be  learned 
from  an  examination  of  Figures  2  0  and  21. 

Figure  20  illustrates  a  point  of  great  opera- 
tional significance o   It  is  seen  that  if  destroyers  were 
designed  so  as  to  tolerate  beam,  quartering  and  following 
seas  operation  from  a  rolling  and  maneuverability  stand- 
point, the  stern  accelerations  could  be  greatly  diminished 
and  thus  the  effectiveness  of  the  DD/VTS  system  is 
greatly  enhanced  in  any  given  sea  state.  A  step  in  the 
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right   direction   can  be    taken   in   the    installation   of 
anti-rolling  fins    on  destroyer-type    ships,    because   as   a 
consequence    of   being   able    to   tolerate   beam  seas    opera- 
tion,   Figure   20  shows   they  could   diminish  stern  accele- 
rations  to   between   50  and    60  percent    of   the   head   seas 
accelerations . 

As   a   slight  digression,    I   see   no  reason  why 
the    same   reasoning  cannot    be   applied   to   problems   associa- 
ted with   bow  movement    such  as   slamming,    shipping   of   green 
water,    degradation   of   bow   sonar   performance   and  general 
ability  to  maintain  high   speed   in  rough  seas. 

Figure   21  illustrates    the   fact   that   ship  length 
has   a  greater   percentage  effect    on   stern  acceleration  in 
small  waves    than   in   large    ones.      Consequently,  increased 
ship  length  pays   significant   dividends  during  the   majority 
of  the   time    because   the    small  waves   are   much  more    prevalent 
than  the    large    ones,    but   pays   lesser  dividends  when  con- 
sidering relative    performance    in   severe   waves. 


66 


UTILIZATION   OF  tJ1"  ACCELERATION  DATA* 

A  crucial  question   is  how  to  most   effectively 
make   use   of  data  such  as  that  plotted   in  Figures  20  and 
21,     By  the   means  discussed   in  Appendix  A,    it   is   possible 
to  predict  from  the  TJE  value   the    probability  of   occurrence 
of  any  particular  value   of  stern  acceleration* 

However,   it   is  highly  desirable   to  simultaneously 
consider  the  normal  distribution  of  "uB  values    (and  thus 
the   Rayleigh  parameter  R)   and   the   Rayleigh  distribution 
of  acceleration  single  amplitudes*     The   theory  for  this 
has   been  developed  by  Nordenstrom  in    (12).     He  assumed 
that    the   parameter  R  in  the   short   term  Rayleigh  dis- 
tribution  is  normally  distributed  and   proceeds  to  obtain 
overall  distributions  by  numerical  integration  of  the 
convolution  of  the  Rayleigh  and  normal  distributions*     A 
sample   of  the   results   of  this  numerical  integration  can 
be  found   in  Figure  22  with  suitable  nomenclature   changes 
made   to  adapt   it  to  the   present  work*     One  can  enter  this 
diagram  with  (Ju/Au  and  for  different  values   of  N   (number 


A  detailed  discussion  of  the  1JE  parameter  is  given  in 
Appendix  A* 
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of  cycles)   obtain  a  factor,   which  when  multiplied  by ''u 
allows    one    to  derive    curves    of  Ax/g   such  as   those    shown 
in  Figure   23   for  the   head   seas   and   570  ft   ship   length 
case  • 

It   must   be    pointed    out   here    that   although  these 
curves    in  Figure   23   could   be    termed    short-term  distribu- 
tion  curves    in   accordance    with   Jasper^   definition    (11), 
I  will  not   do    so,    because    it    is   absurd   to  consider   that 
millions    of   cycles   could   take    place   during  a   short   term 
of   say,    twenty  minutes.      However,    I  do  not   choose    to   call 
them  long-term  distribution   curves,    either,    because    this 
implies   a  natural  vaiation   of   ship   speed   and  heading, 
which   is   not   the    case   here.      These    curves   will  be    con- 
sidered  "long-term"    only   in   the    sense    that  each  curve 
is    the    sum   of  many   short-term  distributions,    keeping    In 
mind    the    assumptions    of   Invariant   ship  speed   and  heading 
that    I  have   been   forced    to  make,    the   former   because    of  a 
paucity   of  model  test   data   and    the    latter   because    of   a 
lack   of   information   available    as    to   probability  distribu- 
tion  of  heading   to   the    sea0      This   matter   Is   discussed 
further   in   the   next   section* 

A   partial  answer  to   this   dilemma   is   to   consider 
the    probability   of   occurrence    of   the    five    sea   states 
shown.      It   is    then   possible   to  weight   these   five   curves 
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by   the   probability  of   occurrence    of  each  sea   state   and 
arrive   at  a   single   curve   that   is   at   least   one    step 
closer  to   being  a  true    long-term  distribution  curve o 

The   following   sea   state   distribution   from   (14) 
for   the   North  Atlantic    (40-50°  North  Latitude)  was  used 
as  an  example  : 

HW3   in  feet  Probability   of   Occurrence 

Calm  1.3% 

less   than  4,92  32.7 

4.92    -   9.84  44.1 

9.84   -   14.77  14.9 

14.77   -   19.69  3.9 

19.69   -  26.20  2.1 

greater  than  26.20  1.0 

l00.($ 

Plotting  this   information   on   probability   paper 

gives    one   a  continuous  distribution  of  HW3.      By  breaking 

this   continuous  distribution  up  into  discrete    increments 

oentered   on  the  H1/3  values  being   used  for   prediction 

purposes    (6,11,17,25   and   37  feet),    it   is   reasonable,   but 

not  mathematically  rigorous   to  state   that: 

Hv/3=6  ft   is  representative   of  what   occurs   70.0%   of  the    time 

*      11  "  22.0  " 

"17  "  5.5  " 

"25  "  2.0  " 

"37  "  .5  " 

166.  <$ 

These    statements  are   not  mathematically   precise 

because   the    ship  response  due   to  a  discrete  S1/3   is   said 

to  be   equivalent   to  the    ship  response   resulting  from  a 

range    of  H3/3  values  above   and  below  that  discrete  value, 
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and   this   is   obviously  affected  by  the   size   of  the   Hw3 
increment   between  discrete   values   of  HWj.      The   larger 
the    increment,   the   larger  would   be   the  deviation  in  "VJB 
values;    larger  than  the  deviations   caused  merely  by 
variation  in  sea   spectral  shape. 

However,  in  the  analysis  undertaken  here,  the 
degree  of  approximation  seems  reasonable  with  general 
regard  ■■  to  accuracy  of  the   predictions. 

The   five   curves   in  Figure  23  were   weighted  at 
various   levels   of  acceleration  by  the   above   probabilities, 
and   the  result  was   the   curve   shown  in  Figure  24  which  I 
choose  to  call  a  quasi-long  term  distribution  ourve   for 
a  particular   length  ship  under  the   assumed  conditions 
of  ship  speed  and  heading.     When  statistical  information 
is  available   regarding  speed  and  heading  distributions, 
inclusion  of   this  information   in  the   analysis  will  make 
a  curve   such  as   in  Figure  24   completely  general,   and 
then  it  could  honestly  bee  called  a  long-term  distribu- 
tion  curve • 

Fortunately,  however,   a  quasi-long  term  distri- 
bution curve   such  as  that   in  Figure  24   is   presently 
useful  for  systems  analysis  and  design  beoause   the   con- 
ditions  of  8 peed  and  heading  used   in   this  work  are   such 
as  to  yield  conservative  results.     In  other  words,  more 
work  and   inoreased  knowledge  will  show  us  how  much 
lower  our  predicticns   can  safely  be.     Here   again  this 
is  given  more  disoussion  in  the   next  section. 
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A  DEMONSTRATION  OP  THE  SYSTEMS  VIEWPOINT 

With  predictions  available  such  as  those 
developed  in  the  preceding  section,  it  is  possible  to 
make  an  evaluation  of  system  effectiveness,  keeping  in 
mind  the  assumptions  as  to  ship  speed  and  heading  that 
were  made.   The  significance  of  these  assumptions  will 
be  discussed  further  later. 

In  the  case  under  study,  the  system  is  the 
DD/VDS  system.  Destroyer  vertical  stern  acceleration 
is  no  doubt  the  most  important  single  parameter  affecting 
operation  of  the  DD/VDS  system,  given  an  otherwise  oper- 
able system.   The  structural  ability  of  the  cable  to 
remain  intact  in  rough  weather  is  a  go/no-go  situation 
dependent  on  vertical  stern  acceleration.   Vertical 
stern  acceleration  also  bears  on  understanding  and 
solving  the  VDS  control  problem  and  on  the  design  of 
the  towing  device. 

A  systems  analysis  has  been  done  herein  using 
vertical  stern  acceleration  as  the  significant  parameter- 
other  parameters  could  be  included  in  the  same  manner, 
if  suchum  deemed  necessary,, 
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There  are  two  basic  areas  of  interest  in  this 
analysis  of  DD/VDS  operational  effectiveness  based  on 
vertical  stern  accelerations;  they  can  best  be  posed  as 
questions : 

(a)  in  the  case  of  an  extant  system,  where  a 
design  assumption  as  to  vertical  stern 
acceleration  has  been  made,  what  should 
the  upper  operational  limit  be,  with 

respect  to  weather,  in  order  to  give 

« 

reasonable   assurance    that   the   cable   will 
not   be    snapped  and  that   the   VDS  might   be 
controllable,    that   is    to   say,    in   order 
that    the   design   assumption    of   vertical 
stern  acceleration  will  not   be   exceeded? 

(b)  in  the   case    of  a  system  under   original 
feasibility   or  design   study,   what   value 
of  vertical  stern  acceleration   should   be 
used  as   a  design  requirement,    so  as   to 
insure    that   the   system  will  be    opera- 
tionally effective   for  an  acceptable   per- 
centage   of  the    time    on  a   long-term  basis? 

For   the   extant   system,    the    curves   in  Figure   23   are 
the   most   useful.      Say,   for  example,    that  Ax  =  1,5  g's  was 
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used   as   a  design   assumption  and    therefore    the   design  was 
based   on   this   figure.      Let   us    say,    further,    that    it    ia 
desired    that    the    probability   of  exceeding    1.5   g   be    .001% 
or   less.      Figure   23   tells   us   that   in    seas   characterized 
by  H^/3   =  25  feet,    the    probability   of  exceeding   1.5  g   ia 
slightly  less   than    . 001#,  which   is   satisfactory.      However, 
in   the    case    of   seas   characterized   by  Hwg   =  37  feet,    the 
probability   of  exceeding   1.5  g   is   about    .095$   ,wnich  ia 
too  great.      Therefore,   we   would    specify  that   this   system 
avoid    operating   in   seas  having  HW3  greater   than  about 
25  feet.      Since    seas   of   this   severity  might   be   expected 
to   occur   only  two   or   three    percent    of  the    time,   this 
system  could  be   expected   to   be    operationally  useful 
97-98   percent    of   the    time,    considering   the    restriction 
Imposed      and    the    probability   of   failure    that  was 
specified   as   being  acceptable.      The    numbers   used   in   this 
example   are   for   illustration   only  -   if  a  smaller   prob- 
ability of  failure   is   desired,   for  example,    a   lower 
upper  sea   severity   limit  would  have   to  be    imposed,    thus 
reducing  the    total  percentage    of  the    time   that   the    system 
could   be    operated.      This   is    obviously  a   situation   involving 
various   trade-offs   and    one   must    pick  the   combination  best 
suited   to  his   needs. 

In  the    case   of  a  new  design,    It   can  be    seen  that 
Figure   24   is   the   most   useful,    if   it   is  desired,    as   it 
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certainly   seems    it  would   be,    to  avoid   any  sea   severity 
upper   limitation    on   the    operation   of    the    system.    Once 
an  assumption    is   made    as   regards    likely   areas    of   opera- 
tion,   a   curve    such  as    that    shown   in   Figure   24    can  be 
used    to  arrive    at   a  reasonable    value    of   stern  accelera- 
tion for   use   as   a  basis  for   design. 

For  example,    let   us    consider  a  DD/VDS   system 
that    is   at    sea   for  40   percent    of  a   ten-year    life    span 
(the    ship  will   last    longer,    but   most,   likely   the    weapons 
system  will   be    changed),    and    assume    that    the    VDS   is 
streamed   for  half   the    time    tile    ship   is   at    sea.      One 
could  expect   about   9  million  cycles    of    stern   accelera- 
tion during   that    period    of   time,    based    on   an     overall 
period    of   stern  motion    of   about    7   seconds.      Entering 
Figure   24   with  N  =  9  million,    one   reads    off   about    1.8  gfs 
as  the    level  of   stern   acceleration  above   which  we   can 
expect   only   one    of   the   9  million   cycles   to   occur.      The 
corresponding   probability   of  failure    if  we    use    1.8  g's 
as   a  basis   for  design   is   about    .00002    percent. 

It   must   be   emphasized   again   at   this    point,  as    in 
the    previous    section,    that   curves    similar   to   those    being 
discussed   are    only  as   good   as    the    assumptions   made    in 
deriving   them.      In   the    specific    case    under   consideration 
in   the    present   work,    the    curves    in   Figures   23   and   24    are 
based   on  a   single    ship   speed    of   17  knots   and,    of   course, 
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the   system  will  not   spend   its  entire    lifetime    operating 
at    17  knots,    although  some    weight   may   be    given   to    the 
fact   that    17  knots   could    be   considered    a   typical   speed 
for    operation   of   the   DD/VDS   system.      Another   fact    that 
must   be    kept    in  mind   regarding   these    particular    figures 
is    that    they   are   for  head    seas    operation  and,    here    again, 
the    system  will  not   spend    its    lifetime    operating   in  head 
seas.      As   a  matter   of   fact,    as   discussed    previously,    it 
would   considerably   lower   the    levels    of   stern  accelera- 
tion  if  we   were    to  design    our  destroyers  with  anti- 
rolling  devices    so  as    to  be    able    to  avoid  head   seas 
operation.      Therefore,    the    curves   given  here   will  yield 
conservative   results    in   any  case. 

My  primary   purpose   in   this    section   is   to  present 
a  method   for   systems   analysis    that   can  be   quite    useful 
in  evaluation   of   system  effectiveness  when   the    input   to 
the   analysis    is   reliable    in  magnitude   and    sufficient    in 
scope.      At    present,    I   consider   the   available    input   in- 
formation  marginal  in   some   respects,    namely: 
(a)      probability  distributions    of  ship 
heading   to   the    sea  a?e   not    presently 
available   and    should    therefore   be 
given   some   attention.      It   seems  evi- 
dent   to  me    on  trie    basis    of  four   years 
at   sea   on   board  destroyers   that   as 
weather   becomes   more    severe,    the 
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tendency  to  steer   so  that   seas   are    on 
the    bow   increases   because    of   a  desire 
for  minimum  motions   and  maximum 
maneuverability.      It   is   due   to  the 
foregoing  that    I  assumed  head   seas   as 
being    slightly  more    conservative    than 
an  assumption    of  bow   seas.      However, 
the   more   widespread   use    of   such  de- 
vices  as   anti-rolling  fins  will  no 
doubt   change   the    situation, 

(b)      it   would  be   highly  desirable    to    include 

the    effect   of   speed    in   these    predictions 

even  though  it   probably  has   a   second   order 

effect    compared   with   ship   length  and    sea 
« 

severity.      At    present   there    is   no   model 

test    series    available    that    covers   all 
headings    and   speeds.      Perhaps    analytical 
derivation   of  the   RAOfs,    such  as   is   now 
being    investigated   at   Webb,   would    solve 
this   dilemma,    by  replacing   the   model 
testing   process   with  an   analytical   pro- 
cedure   than   can   be    done   with  a   computer. 
One   further  point   that  must   be   discussed   is   that 
of  attempting   a   complete    systems   analysis.      There   are 
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many  factors  that  affect  a  system  as  complex  as  a  ship 
and  its  equipment,  and  I  have  examined  only  one  sub- 
system of  the  larger  overall  system.   The  reliability 
of  the  machinery  in  its  task  of  motivating  the  ship 
during  VDS  operations  at  the  desired  speed  certainly 
affects  the  effectiveness  of  the  DD/VDS  system,  but  the 
inclusion  of  factors  such  as  this  is  left  to  the  opera- 
tions analyst  charged  with  maintaining  the  macroscopic 
view,  more  commonly  referred  to  an  the  "big  picture." 
It  is  my  hope  that  the  information  developed  in  this 
thesis  regarding  destroyer  stern  acceleration  might  aid 
him  in  his  task. 
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CORRELATION   OP    PREDICTIONS   WITH   FULL-SCALE   MEASUREMENTS 

Aa   a  result    of   interest    in  VDS   systems,    the   U.S. 
Navy  Underwater  Sound   Laboratory   in  New   London,    Connecti- 
cut,  has  recently  embarked   on  a  program  of   collecting  a 
statistical  sampling  of  vertical  stern  accelerations    on 
destroyers   at   sea    (17).      For  this   purpose,   G-iannini 
Controls   Corporation  Model  2452   statistical  accelero- 
meters  have    been   installed   on  various   destroyers   at   the 
stern,   as  well  as   other   locations. 

Up   to   the    present   time,    these    accelerometers  have 
been   installed   so  as   to   detect   instantaneous   values   of 
positive   acceleration   that   exceed   the   following  four 
thresholds:      1.20g,    1.35g,    1.50g  and   1.65g,    and  record 
on   the   four   counters,    the  number   of  exceedences  above 
each   of  the    thresholds© 

It  must   be    pointed   out   here    that   since   these 
accelerometers   are   designed   to  be   used  horizontally  as 
well  as   vertically,   the    zero  reference   for  measuring 
vertical  accelerations   is   one    g  rather   than   zero  g,    the 
latter  being   used   throughout   this   thesis.      As  a  conse- 
quence,   the    above    thresholds   become    .20g,    .35g,    .50g 
and    .65g   in   the    zero  g  reference    system  for  vertical 
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ace- lerations .   No  attempt  will  be  made  here  to  debate 
the  pros  and  cons  of  the  two  systems  -  suffice  it  to  say 
that  there  are  arguments  for  both. 

As  of  the  present,  there  has  been  collected  a 
small  amount  of  this  statistical  data  (and  to  my  know- 
ledge it  is  the  only  such  data  available)  and  it  is  of 
considerable  interest  to  compare  it  with  the  predictions 
made  in  this  thesis. 

The  data  received  from  the  Underwater  Sound 

Laboratory  is  as  follows: 

USS  MOAhB  (DP 6 95) 

Number  of  accelerations  above  .20g  =  708 

.35g  =  239 
.50g  =  44 
.65g  =   4 

Total   time    of   test    =  226  hours 
Estimated  number   of  cycles   for  6   second 

period   of  motion  =  135,600  cycles 
Weather      l/6   -  6  foot  wave   heights    observed 
Headings   and   speeds   -  various 

The    probability   of  exceeding  various    levels    of 

acceleration,    therefore,   would   be   as   follows: 

Probability   of  exceeding    .20g  =  708/135,600  =  .530% 

.35g  -  239/135,600  *  .176 

.50g  ■     44/135,600  ■  .0325 

.65g  =        4/135,600  =  .00294 


USS    PUTNAM    (DP    757) 

Number  of  accelerations  over   .18g  =  6009 

,38g  =  1529 
.46g  =   149 

.67g  ■   17 
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Total  time   of   test   612  hours 

Estimated  number   of  cycles  for   6  seoond 

period  of  motiaa  ■  367,000  cycles 
Weather  1-10  foot  waveheights  observed 
Headings   and   speeds   -  various 

Probability  of  exceeding   .18g  =  1,64% 

.38g   m     .43 
,46g   =      .041 
,67g  =      .0046 

USS  H.    J.   ELLISON    (DP   864) 

Number   of  accelerations   above    .19g  =  1229 

.34g  ■     156 

,48g  -        18 

•  66g  =  1 

Total  time    of  test      629  hours 
Estimated  number   of  cycles   for   6   second 

period   of  motion  =  377,000  cycles 
Weather  0-12  foot  wave  heights   observed 
Headings   and   speeds   -  various 

Probability  of  exceeding    .19g   =.34% 

.34g  =.043 
.48g   =.0049 
.66g   =.0003 

The   Moale   and   Putnam  are   destroyers   of  the   Sumner 
Class  whose    length  on  the   waterline    is  369  feet.     The 
Ellison  is  a  Gearing  Class  Destroyer  whose   length   on  the 
waterline   is  383  feet. 

The    probabilities   of  exceeding  the   four   levels 
of  acceleration  for  each  ship  were    plotted   in  Figure   23  (Pg.70) 
to  compare   them  with  predicted  values. 

Since   the    predicted  curves   are  for  head   seas 
and   17  knots   ship  speed,    one   would   expect  them  to   be 
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high  relative   to  the    observed   data  which   is   for   various 
headings   and    speeds.      Therefore,    since    the    observed   data 
was   taken   in   seas   in  the   neighborhood   of   0-11  feet   in 
height,    one   would  expect  the    observed  data  to  plot   some- 
where   below  the   Hw3   =  11  feet   curve.      As   can  be    seen   in 
Figure   23,    this    is   exactly  where   the    observed   points    plot. 
The   dashed    line    shows    the    trend   for  the    two  3  69  foot 
ships   and    the   dot-dash  line    shows   the    trend   for   the   383 
foot   ship.      The    longer   ship  exhibits   a   tendency  toward 
lower   stern  accelerations  which  agrees  with   the   trends 
of  stern  acceleration  with   ship   length  previously  dis- 
cussed. 

The   fact  that   the    lines  drawn   through  the    ob- 
served  data  have   a   larger   slope    than   the    individual  H3V3 
curves   seems   to  reflect    the   fact   that   the    observed  data 
included   the   effect    of  a  range    of  weather  conditions. 
The    curve    of  predicted  values   in   Figure   24  has   a   larger 
slope   and    it   likewise   includes   the   effect   of  a  range    of 
weather  conditions. 

This  section  is  intended  in  no  way  to  be  a  com- 
prehensive correlation  of  observed  data  with  predicted 
data.  Indeed,  much  more  observed  data  is  needed  so  one 
might  cull  it  for  information  as  to  accelerations  for  a 
particular  heading  and  speed.  However,  the  information 
presented  here    tends   to   increase   one's   confidence   in   the 
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quantitative    and   qualitative    aspects    of   the    stern 
accelerations    predicted   by  the    method    presented   in   this 
thesis. 
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CONCLUSIONS 

1*     The   method  presented   in  this   thesis   pro* 
vides  the  system  designer  with  a  powerful  technique   — 
where    ship  response   to  the   sea  is  a  pertinent  factor  — 
for  analyzing  on  paper  the   worth  of  a  proposed  system  while 
still  in  the  conceptual  stage,' 

2.  The  available  full-scale   stern  acceleration 
data  tends  to  confirm  the  qualitative  and  quantitative 
aspects   of  the   predictions  made   by  the   method   presented 
herein. 

3.  Prom  the   standpoint   of  diminishing  stern 
acceleration,    ship  heading  to  the    seas   is  the  most    impor- 
tant  single  controllable  factor  and,   therefore,  destroyer- 
type  vessels  with  installed  equipment   sensitive  to  stern 
accelerations  should  be  designed  so  as  to  be  able   to 
tolerate  rough  seas  operation   in  beam,  quartering,  and 
following  seas*      Installation  of  anti-rolling  devices   Is 

a  step  in  the   right  direction  by  enabling  the   ship  to 
operate   in  beam  seas  and  thereby  paying  a  dividend  of  50- 
60  peroent  diminution  in  stern  acceleration* 

4.  Ship  length  has  a  greater  percentage  effect 
on  stern  acceleration  in  light   seas  than  it  does   in 
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heavy  seas.      Consequently,    increased   ship    length  pays 
significant   dividends  during   the  majority  of   the  time 
since   the    light   seas   are  much  more   prevalent   than  are   the 
heavy  ones*     The  effect   of   ship   length  is  not  so  marked 
when  considering  relative   performance    on  a  percentage 
basis   in  severe   seas* 

5*     Sea  spectrum  shape  has  a  significant  effect 
on  stern  acceleration  response   and  should  therefore   be 
accounted  for  in  the   method  used  to  predict  that  response. 
The   method  used  here   seems  to  be  effective    in  that  regard* 

6*     As  demonstrated   in  this   thesis,   the   system 
approach  seems   particularly  useful   in  certain  aspects   of 
ship  design,    since  a  ship-related   system  may  most    logically 
be  ccnsidered   in   terms   of   its   intended  mission,  for  example* 
the    location  of  enemy  submarines* 

7*      In  a  problem  such  as   the    one   studied  here* 
the   use   of  a  high  speed  electronic   computer  virtually 
frees  the  designer  from  worries  regarding  time   required 
for  calculations  and  allows  him  to  spend  full  time   in- 
suring the   soundness  of  the  analytical  processes  being 
used* 

8*     The  generality  of  the   analysis  developed 
here    insures  that  a  similar  method  can   be   utilized  to 
investigate   other  system  characteristics  related  to  sea- 
keeping*   such  as  shipping  water  over  the   bow,   sonar  de- 
gradation and   slamming* 


87 


9,      The    analysis    used    in   this    thesis    shows    that 
the    behavior  and   effectiveness    of   a   ship-related   system 
is    ideally  suited   for  expression   in   probabilistic    terras. 
However,   both   ship  designers    and    operators   will  need   to 
develop  a  whole    new  framework   of   thought   and    reference 
as   to   acceptable    probability   levels    of    system  failure 
and/or   loss.      Indeed    it    puts   a   designer  face    to  face 
with  his    responsibilities   when  he   has    to  express    loss    of 
resources,   human  and    otherwise,    in  mathematical  terms 
that    are    uncomfortably   precise. 
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RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  Investigate  means  that  would  allow  destroyer- 
type  vessels  to  tolerate  operation  in  rough  quartering 
and  following  seas.   The  obstacles  to  be  overcome  include 
considerations  of  stability  and  maneuverability. 

2.  Continue  work  presently  being  done  in  the 
development  of  theory  for  analytical  derivation  of  res- 
ponse amplitude  operators,  including  experimental  verifi- 
cation. 

3.  Investigate  probability  distributions  of  ship 
heading  to  the  sea  under  various  conditions. 

4.  When  data  becomes  available,  undertake  a 
comprehensive  correlation  between  full  scale  stern 
acceleration  measurements  and  the  predictions  contained 
herein. 

5.  Give  increased  emphasis  to  collecting  more 
seaway  measurements  useful  for  ship  design,  preferably 
directional  data  measured  with  a  buoy  rather  than  a 
shipborne  recorder.   In  this  regard,  effective  liaison 
between  ship  designers  and  oceanographers  is  of  the  ut- 
most importance. 

6.  Investigate  stern  vertical  velocities  and 
displacements  under  various  conditicns  of  weather,  ship 
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length  and   heading   using   a   method    similar  to   the    one 
presented  here.      Also   investigate    the   effect    of   ship 
speed    on    stern  vertical  acceleration,    velocity   and   dis- 
placement . 

7.      Investigate    the    application,    to  the    VTS 
cable    stress    situation,    of  an    Instrument    similar  to   the 
flexing   stress   monitor   described    in    (18).      This    instru- 
ment   could    be    used    to   compute    and    present   maximum   probable 
cable    stress   during   the   ensuing    15  minutes    or    so   of 
operation   and    thus   give    the    ship's    personnel  an   objec- 
tive   basis   for  making   speed   and/or  heading   changes   when 
the    probable    cable    stress    becomes   excessive. 
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APPENDIX  A 


DISCUSSION   OP  THE  IS   PARAMETER 


A  record    of  any  wave-caused    phenomenon!  in 
irregular   seas   such  as   stern  acceleration  would   look 
something  like 


Snp* 


where  x>      is   stern  acceleration  measured  from  the   mean 
value   at  equal  or  random  time   intervals   and  yA    is   peak  to 
peak  stern  acceleration  measured  wherever   it    occurs* 

It   is  now  widely  accepted   in   the    literature*  that 
the  xx    *»   are   normally  distributed   and   that   the   yx's   are 
Rayleigh  distributed   on  a  short-term  basis,     A  short- 
term  distribution  if  applicable    so  long   as   uniform  condi- 
tions  of   sea,    ship  speed   and   course    persist.      This    is   to 


*See  for  example,    (11)  and   (13). 
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be  distinguished  from  a  long-term  distribution  which  allows 

the  sea,  ship  course  and  speed  to  vary  over  a  range  of 

conditions.  Jasper  introduced  these  definitions  in  (11) 

and  they  are  used  by  many  investigators,  including  Nordenstrom  (12). 


(Peak  to  peak) 


Rayleigh 


:  f  (y)  -  <* 


-y2^ 


(Sucessive 

Ordinates ) 


Normal: 


(TPW 


It  is  seen,  then,  that  the  parameters  necessary  to  define  the 

2 
above  distributions  for  a  particular  case  are  R  and  <T~  •  Fortunately, 

it  is  possible  to  get  these  parameters  from  the  Integrated  Response 

Spectrum  which,  through  the  superposition  concept,  has  certain 

mathematical  properties  in  common  with  the  irregular  record  shown. 


»  for  brevity,  the  normal  probability  distribution 
function  is  often  referred  to  as  N(//,j^)c 
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It  has   been  shown  in   the    literature4**  and   is   now 
widely  accepted   that,   for  the    type   of   spectrum  being   used 
in  this  work, 

R  a  8  6*2  -  4B  (13) 

Thus   it   la   seen   that  R  and  6" 2  are   known   once  E  is  known 
from  the    Integrated  Response  Spectrum, 

It   can  also  be   demonstrated,   using  the  concept 
of  maximum  likelihood  estimation,   that 

R    »>      7k2 

<m    n 
and   since  root  mean  square    (of  peak  to  peaka )   is  defined 
as 


to) 

7 


RMS     f'peakry      *  ]j   L.y^2 
one   may  relate  RMS(ptp)   to  R  and   thus  E  as   follows 

RMS(ptp) 


=  lLLjk£      =t/r  -  2lfe 


RMS  of  peak  to   peaks   is   commonly  referred  to  as  RMS  of 
double   amplitudes.     Therefore  RMS  of  single   amplitudes 

RMS/single     \   *  RMS    .    /2   =  7JI" 
(amplitude^  pT;p 

Thus   it    is   seen  that   the  ifE  values   plotted   in 

the   various  figures   can  be    thought    of  as  being  the  root 

*S-   (13) 
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mean   square    of   single   amplitudes    of   stern  acceleration 
and   such  is  the   case   throughout   this    thesis. 

All   of   the   above   relationships   were    originally- 
developed   using   an   irregular  record    of  wave   height   and 
the    corresponding    sea   spectrum,    but    I  have   discussed 
the    situation   in   terms    of   stern  acceleration  because 
ship  characteristics   dependent   upon  waves   such  as    pitch 
and  heave   have    been  found   to  behave    in  the    same   fashion 
as   the   waves. 

Since   most    often   one    is   most   interested   in   single 
or  double   amplitudes    of  motion,    it    is    common   practice 
to   use    the    parameter  E  which  defines   the   parameter  K 
and   thus    the  Rayleigh  distribution   of  amplitudes    to 
arrive    at,    for  example,    the   average   highest   expected 
value    of   stern  acceleration   single    amplitude    in  20  cycles 
from 

y  Ax/gh/gQ  =  1.87  t/b 

Multipliers   for   other  numbers   of  cycles  N   (where 
N>50)   may   be    computed   to  within  3%  accuracy  using   the 
following  formula  from   (10): 

(In  J)  +  y  (In  N)"1^      where  r=  Eulers  Constant   = 

2  .5772 

Examples   are   as  follows: 

N=n umber   of   cycles  Average    of  highest  expected 
m         values    of   single   amplitude 

50  2.12  7/1 

100  2.28  TfE 

200  2.43  iff 

500  2.61  tfS" 

1000  2.74  ifE 
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For   the   reader's    convenience,    the   following 
are    also  quoted:  (6) 

Most   frequent    single    amplitude  =  0.707  liE 
Average    of  all   single    ampli- 
tudes =   0.886  "tfl 
Average    of   l/3  highest   single 
■                                           amplitudes  =  1.416  Z)B 
Average    of   l/lO  highest    single 

amplitudes  =  1.800  ZlB 
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APPENDIX  B 

A   UNIFORM  DEFINITION    OF   SEA   STATES 
For  about  fifteen  years,   a  world-wide   sea  state 
scale    (World   Meteorological  Organization  Code   75)    (15) 
has   been   in  general  maritime    use   by  both  naval  and 
merchant   ships  for  the   collection  of  weather  data  at 
sea* 

It  seems   logical,    therefore,    in  the   interest   of 
meaningful  full-scale   correlation,   that  all  our  engineering 
predictions  for  ultimate  application  to  the   ships  at   sea 
should  be    in  terms   of  WMO  Code   75,  which  is  quoted  below 
for  the   reader's  convenience: 

STATE   OF  SEA -WIND   WAVES    (WMO  COIE  75) 
Code  Description  Height    in  Feet  # 

0  Calm-glassy  0 

1  Calm-rippled  0-1/3 

2  Smooth-wavelets  l/3-l  2/3 

3  Slight  1  2/3-4 

4  Moderate  4-8 

5  Rough  8-13 

6  Very  Rough  13-20 

7  High  20-30 

8  Very  High  30-45 

9  Phenomenal  over  45 


*It   is   proposed  herein  to  consider  this  column  as   being 

Vs. 
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In  spite   of  the   above,   other  sea   state   scales, 
such  as   the   one  reproduced   in  Figure   25,   are    in  wide   use 
by  designers,   presumably  because   of   the  mathematically 
convenient   use   of  significant  wave  height    (H1/3)  versus 
sea  state  designations*     However,    one   simple    and  defensible 
assumption  would  remove   the  need  to  use  anything  but  WHO 
Code   75.  This   is  discussed   in   the   two  following  paragraphs. 

Both  references    (6)  and   (16)  agree  with  what  many 
oeeanographers  seem  to  feel  in  that  at -sea  wave  height 
observers  tend  to  overestimate  wave  height*     Reference 
(16)  even  goes   so  far  as  to  instruct  the   observer  to 
record  the  average  of  the  estimated  heights  of  the 
larger,   we  11 -formed  waves* 

Consequently,   it  seems  to  make  a  great  deal  of 
sense   to  specify  the  assumption  that  observed  wave  height 
equals   significant  wave  height.   H^,  where  fi^  is   the 
average  height  of  the   1/3  highest  waves   in  a  seaway*     This 
would  seem  to  bridge  the   gap  between  the  designers  who 
have  been  prone  to  use   such  scales  as  the   one  shown  in  Figure 
25   because   of   its  mathematically  desirable   use   of  R1/5  and 
the   ships  at  sea  who  have   been  making  observations  based 
on  WMO  Code  75* 

It  might  be   pointed  out  that   the  Hw,  values   of 
6,    11,   17,  26  and  37  feet  used   in  this   thesis  for  pre- 
dictions of  stern  acceleration  were  chosen  as   being 
representative   of  WMO  Code   75  Sea  States  4,   5,   6,   7  and 
8  respectively* 
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